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CROSS-REFERENCE TO RELATED APPLICATIONS 

This application is a continuation-in-part of U,S. Application Serial Number 10/389,183, filed 
March 14, 2003. 

FIELD OF THE INVENTION 

One embodiment of the present invention relates to a method for constructing a circuit for 
controlling an electromagnetic actuator. 

Another embodiment of the present invention relates to a method for designing a circuit for 
controlling an electromagnetic actuator. 

For the purposes of the present application the term "physically remote" (e.g., in the context 
of a coil being physically remote from an electromagnetic actuator) is intended to refer to the fact 
that the electromagnetic actuator and the coil may be electrically connected but that any direct 
magnetic interaction between the two is negligible. 

Further, for the purposes of the present application the term "theoretical" (e.g., in the 
context of a theoretical coil) is intended to refer to the fact that the theoretical coil does not exist in 
the physical sense. 

BACKGROUND OF THE INVENTION 

In general, a solenoid converts electric energy into magnetic flux, release of which is 
transferred into linear mechanical motion of a plunger installed in the center of a C-frame solenoid, 
a D-frame solenoid, or a tubular solenoid (as shown respectively in Fig. 1 A, Fig IB, and IC). 
Current flow / through the solenoid coil winding with inductance L creates magnetic energy 

£■ = which produces an attraction force Fmag between a movable plunger and a fixed stop. 

Solenoids typically have a working, or variable, air gap between the plunger and the stop, as well as 
a fixed air gap between the outside diameter of the plunger and either its frame or mounting 
bushing. To complete the magnetic circuit, the magnetic flux lines flow through either air or the 
metallic frame through the stop, the plunger, the frame or the mounting busing of a tubular solenoid 
and return to their point of origination. 

The performance of a solenoid is dependent on numerous parameters, including, but not 
necessarily limited to, its physical size, the wattage 

applied, duty cycle, ambient temperature, its coil temperature due to heat rise, the coil 
ampere-turns (NI where / and N are current and coil turns respectively), solenoid orientation, cross 
sectional area of the plunger, the coil winding and the plunger and stop geometry. Fig. 2 illustrates 
typical force-stroke relationships for different geometries of plunger and mating stops of a D.C. 
solenoid. 

Typically, the greater the holding force of a given plunger and stop geometry, the lesser the 
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pulling/pushing force at an extended stroke position. In this regard, the minimum pull/push force 
generated is typically at the extended stroke end where the plunger assembly begins it's lifting 
towards the stop. As the plunger approaches the stop position, the pulling/pushing force developed 
typically increases dramatically, and the slope of the force-stroke curve rises sharply. The 
differential equations for an electrical circuit and Maxwell's equations for dynamics, which define 
the forces according to the current and position, describe the full dynamic or switching response of 
an electromechanical actuator. In fact, there is a certain transient time needed to develop magnetic 
flux and transfer it's energy to mechanical momentum. 

In many applications this intrinsic transient phenomenon may ultimately effects the 
dynamics of other mechanical parts dependent on the plunger position and it's speed. One 
of these applications is related to high-pressure fuel injectors used in direct injection gasoline and 
diesel engines. In internal combustion engines (especially diesel engines) the transient phases, 
including injection, ignition (or auto-ignition) and combustion, have ultra-short time fractions from 
a few tens to a few hundreds of a nanosecond. In this regard. Fig. 3 shows data regarding normal 
heptane reactions starting at 900K and 83 bar in connection with a two stage CI (diesel) combustion 
process. More particularly, Fig. 3 relates to: (a) a first stage including premixed flame (0.03 ms) 
having various short-lived species such as C7 radicals, aldehydes (PAH), and hydrogen peroxide; 
and (b) a second stage including rapid oxidation (0.06 ms) having hydrogen, water, carbon dioxide, 
carbon monoxide, methane, soot precursors, C3-compounds, and C4-compounds. 

Further, Fig. 4 depicts certain ideally targeted or aimed or purposed injection events (e.g., 
hampered by unstably controlled injection shot duration and dwell interval) and Fig. 5 depicts a 
diesel difftision flame in connection with a conventional single long shot per cylinder injection 
(with limited access of air resulting in incomplete combustion). 

Further still, one conventional electronically controlled diesel fuel injector is called an 
"accumulator" type. In these injectors, a nozzle includes an accumulator chamber that is charged 
with fuel under high pressure, which communicates with a nozzle port. An actuating device is 
associated with the injection valve and is moveable within a control chamber that is also 
pressurized with fuel. A valve is associated with the control chamber and is opened so as to reduce 
the pressure and cause the pressure in the accumulation chamber to unseat the injection valve and 
initiate fuel injection. Typically, a main electromagnetic assembly that is contained within the 
housing of the fuel injection nozzle operates the valve. 

Figs. 6A-6D depict four strokes of unit injector ("UI") and unit pump ("UP") operation 
stages. The function of these single-cylinder injection-pump systems can be subdivided into four 
operation stages (corresponding, respectively, to each of Figs. 6A-6D): 

a) Suction stroke . The follower spring (3) forces the pump plunger (2) upwards. The fuel in 
the fuel supply's low-pressure stage is permanently under pressure and flows from the low- 
pressure stage into the solenoid-valve chamber (6) via the bores in the engine block and the 
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inlet (or feed) passage (7). 



b) Initial stroke . The actuating cam (1) continues to rotate and forces the pump plunger (2) 
downwards. The solenoid valve is open so that the pump plunger (2) can force the fuel 
through the fiiel-retum passage (8) into the fuel supply's low-pressure stage. 

c) Delivery and injection stroke (or Prestroke) . An electronically timed signal from the engine 
electronic control unit ("ECU") energizes the solenoid-valve coil (9) to pull the solenoid 
valve needle (5) towards the solenoid valve seat/stop (10). The connection between the 
high-pressure chamber (4) and the low-pressure stage is closed. Further movement of the 
pump plunger (2) causes increased fuel pressure in the high-pressure chamber (4); the fuel 
is also pressurized in the nozzle-needle (or nozzle assembly)(l 1). Upon reaching the nozzle 
needle opening pressure (typically over 300 bar), the nozzle needle (1 1) is lifted from its 
seat and fiiel is injected into the engine combustion chamber. Due to the pump plunger's 
high delivery rate, the pressure continues to increase throughout the whole of the injection 
process (typically up to maximum peak of 1800-2000 bar). 

d) Residual stroke . As soon as the solenoid-valve coil (9) is switched off, the solenoid valve 
(or solenoid valve needle) (5) opens after a short delay and opens the connection between 
the high-pressure chamber and the low-pressure stage. 

Figs. 7A-7D relate to the above-mentioned operating stages of Figs. 6A-6D and show, 
respectively, coil current (Is), solenoid-valve needle stroke {huX injection pressure (pe), and nozzle- 
needle stroke (hn). 

Fig. 8 depicts a wave form diagram associated with operation of a fuel injector nozzle (an 
"accumulator" type injector) under use of two actuating solenoids installed into injector. 

Finally, a number of conventional techniques and apparatuses achieve multiple injection, 
for instance, by using a piezoelectric actuator during individual injection phases or a rapid 
switching on/off of injection events strategy via the electronic control unit. Specifically with 
reference to application of rapidly operating electromagnetic actuators, studies have been carried 
out on variable valve actuators for valve train parts, rather than for high-pressure fuel injectors. 
Related documents include: 1) Robert Bosch GmbH (1999). Diesel-engine management. SAE, 2"^ 
edition, 306 p.; 2) B. Riccardo, C.R.F. Societa' Consottile per Azioni (2000). Method of controlling 
combustion of a direct-injection diesel engine by performing multiple injections. European patent 
EP 1 035 314 A2; 3) N. Rodrigues-Amaya, et. al. (2002) Method for injection fuel with multiple 
triggering of a control valve. Robert Bosch GmbH, US patent 2002/0083919 Al; 4) M. Brian, 
Caterpillar Inc. (2002). Method and apparatus for delivering multiple fiiel injection to the cylinder 
of an engine wherein the pilot fuel injection occurs during the intake stroke. Intentional patent WO 
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02/06652 A2; 5) K. Yoshizawa, et. al, Nissan Motor Co., Ltd (2001). Enhanced multiple injection 
for auto-ignition in internal combustion engines. Us patent US 2001/0056322 Al; 6) Y. Wang et. 
al., Ford Motor Company and K.S. Peterson et. al.. University of Michigan (2002). Modeling and 
control of electromechanical valve actuator. SAE Intenational,SP-1692, 2002-01-1 106, 43-52; and 
7) V. Giglio et. al. (2002). Analysis of advantages and of problems of electromechanical valve 
actuators. SP-1692, 2002-01-1106, 31-42. 

BRIEF DESCRIPTIQN OF THE DRAWINGS 

Figs. 1 A-IC depict, respectively, typical cross-sections (with magnetic flux line patterns) of 
a C-frame solenoid, a D-frame solenoid, and a tubular solenoid; 

Fig. 2 depicts typical force-stroke relationships (curves) for various conical, flat face, and 
stepped conical plunger- stop configurations for a D.C. solenoid; 

Fig. 3 depicts data regarding certain heptane reactions in connection with a two stage CI 
(diesel) combustion process; 

Fig. 4 depicts certain conventional injection events; 

Fig. 5 depicts a typical diesel diffusion flame in connection with a conventional single long 
shot per cylinder injection (with limited access of air resulting in incomplete combustion); 

Figs. 6A-6D depict four strokes of unit injector ("UI") and unit pump ("UP") operation 

stages; 

Figs. 7A-7 relate to each of the stages of Figs. 4A-4D and depict, respectively, coil current 
(Is), solenoid-valve needle stroke (hM), injection pressure (pe), and nozzle-needle stroke (hn). 

Fig. 8 depicts a wave form diagram associated with operation of a fiiel injector nozzle 
example (an "accumulator" type injector) under use of two actuating solenoids installed into 
injector; 

Fig. 9 depicts applied forces at the start and the end of injection according to an 
embodiment of the present invention. 

Fig. 10 depicts a graphic of one example of an I-Function (i.e., If,{t) and its first order 

derivative ^^^^^/^^) according to an embodiment of the present invention directed to a single 

injection event; 

Fig. 1 1 A depicts one example of a secondary coil incorporated into an electric control 
circuit according to an embodiment of the present invention and Fig. 1 IB depicts two associated 
timing scenarios according to an embodiment of the present invention (wherein the top diagram in 
Fig. 11 B indicates charging of a secondary coil simultaneously with injector firing (simultaneous 
charge) and the bottom diagram in Fig. 1 IB shows charging of the secondary before the injector 
firing (pre-charge)); 

Fig. 12A depicts one example of waveform time series for a simultaneous charged 
secondary coil according to an embodiment of the present invention (wherein the bold solid line is 
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a triggering signal controlling injection duration by T2 of Fig. 1 1 A (CD cycle of Fig. 1 IB) and the 
regular solid line is output voltage measured from primary coil) and Fig. 12B depicts one example 
of waveform time series for a pre-charged secondary coil according to an embodiment of the 
present invention (wherein the bold solid line is a triggering signal controlling injection duration 
by T2 of Fig. 1 lA (CD cycle of Fig. I IB) and the regular solid line is output voltage measured 
from primary coil). 

Fig. 13 depicts stable multiple ultra-short injection according to an embodiment of the 
present invention; 

Fig. 14 depicts one example test system configuration used for verification of time response 
dynamics according to an embodiment of the present invention; 

Fig. 15 depicts one example injection system test cell according to an embodiment of the 
present invention, which test cell is used to verify reaction of a fiiel injector connected in series 
with a charged secondary coil (instantaneous fuel flow rate measurements using laser Doppler 
anemometer indicate real fuel dynamics while injection oscillatory flow in capillary quartz pipe). 

Figs. 16A and 16B depict example plots according to an embodiment of the present 
invention of a comparison of different secondary coil ("SC") charging scenarios at the same 
injection condition (Fig. 16A relates to instantaneous volumetric flow rate and Fig. 16B relates to 
integrated injection mass) (flow measurement results); 

Fig. 17 depicts a series of example plots according to an embodiment of the present 
invention of instantaneous volumetric flow rate (top row) and integrated mass (bottom row) time 
series obtained for different charging schemes (i.e., simultaneous charge - 1st column, pre-charge - 
2nd column; and shifted charge - 3rd column) (flow measurement results); 

Fig. 18 depicts one example of controllable high-pressure multiple injection according to an 
embodiment of the present invention; 

Fig. 19 depicts certain injection events associated with one example of an embodiment 
according to the present invention (wherein the injection events are identified with reference to 
certain combustion effects and engine runs/injection strategies); 

Fig. 20 depicts information relating an embodiment of the present invention ~ that is, 
information relating to RL measured (left, primary) and calculated data (right, secondary); 
Inductance and resistance data measured "out of circuit"; L/C meter IIB; L_stray = 2.139 w H, 
R_stray = 0.2-0.3 W; 

Fig. 21 depicts one example of an I-Function arbitrary current trace normalized to unit and 
its first derivative according to an embodiment of the present invention; 

Fig. 22 depicts one example of an I-Function current fitted to certain library rise and fall 
exponential functions according to an embodiment of the present invention; 

Fig. 23 depicts data relating to one example secondary coil driver code (e.g., relating to the 
calculation of certain parameters) according to an embodiment of the present invention; 

Fig. 24 depicts data relating to construction of a current waveform for multiple injection 
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(e.g., associated with an HP Agilent 3481 1 A/331 20A configuration) according to an embodiment 
of the present invention; 

Fig. 25 depicts certain example signals constructed as arbitrary waveshapes (wherein the 
left plot is associated with an original Bosch CRIS injector signal and the right plot is associated 
with a two shot injection signal according to an embodiment of the present invention); 

Fig. 26 depicts one example controllable multiple injection system (applied to a Bosch 
common rail system) according to an embodiment of the present invention; 

Fig. 27 depicts one example measurement setup to verify high pressure multiple injection 
according to an embodiment of the present invention; 

Figs. 28-45 depict the performance evaluation of a multi-burst rapidly operating secondary 
actuator according to an embodiment of the present invention as applied to a diesel injection system 
(of note, this rapidly operating secondary actuator according to an embodiment of the present 
invention may hereinafter sometimes be referred to as "ROSA"); and 

Figs. 46-70 depict the quantification of instantaneous diesel flow rates in flow generated by 
a stable and controllable multiple injection system (i.e., "ROSA") according to an embodiment of 
the present invention. 

Among those benefits and improvements that have been disclosed, other objects and 
advantages of this invention will become apparent from the following description taken in 
conjunction with the accompanying figures. The figures constitute a part of this specification and 
include illustrative embodiments of the present invention and illustrate various objects and features 
thereof 

DETAILED DESCRIPTION OF THE INVENTION 

As required, detailed embodiments of the present invention are disclosed herein; 
however, it is to be understood that the disclosed embodiments are merely illustrative of the 
invention that may be embodied in various forms. In addition, each of the examples given in 
connection with the various embodiments of the invention are intended to be illustrative, and not 
restrictive. The figures are not necessarily to scale, some features may be exaggerated to show 
details of particular components. Therefore, specific structural and functional details disclosed 
herein are not to be interpreted as limiting, but merely as a basis for the claims and as a 
representative basis for teaching one skilled in the art to variously employ the present invention. 

In summary, various embodiments of the present invention relate to electromagnetic 
actuators used to control fuel injectors in internal combustion engines, linear solenoids, and other 
electromagnetic devices (e.g., which convert electric energy into a linear mechanical motion to 
move an external load a specified distance). More specifically, various embodiments of the present 
invention describe the theory, electric circuit, charge time computing code, and engineering 
applications of a secondary coil ("SC") that generates what is herein referred to as an "I-Function" 
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to be used for energizing a first main coil (e.g., installed in a device such as an internal combustion 
engine's fuel injectors). Of note, effects produced by the SC according to the present invention may 
be realized via means taking at least three different forms: (a) an extra, secondary coil installed 
physically remote from the first one (e.g. medium and heavy load solenoids for gasoline and diesel 
engines, for example); (2) an electronic current simulation circuit (e.g. lower load devices, for 
example); and/or (3) a digital/binary code generating an I-Function applied to a desired application 
(e.g., a fuel injector). 

Of further note, three basic problems of mechanic dynamics, induction dynamics, and a 
rapidly operating control unit using an SC are addressed in connection with suppression of any 
transient inertia (delays). In one embodiment the analytical solution is based on a series of 
differential equations. A two-coil configuration of an embodiment of the present invention, for 
example, does not rely upon the physical placement of the second solenoid relative to the first 
solenoid in order to improve valve-lifting response based on the magnetic flux interference between 
the primary and secondary coils. Rather, the present technique realizes an "I-Function" current to 
be applied onto the primary coil. The current may be generated in a secondary coil (which need not 
be physically present in vicinity of the first coil). The secondary coil may be a remote unit that may 
be located away from the first one. The secondary coil may alternatively be presented by a code of 
I-Function induction current to be transmitted and applied. Thus, essentially any desired kind of 
switch on/off process now may be released very rapidly without substantial time lag sensitive to the 
process (e.g., in connection with a combustion process in diesel engines). 

Further, the present invention provides an embodiment in which an electric circuit is 
provided (as well as the code to compute the charging (energizing) time of the SC). In one example 
(which example is intended to be illustrative and not restrictive), the present invention may permit 
injection in a diesel engine in series of pilot and multi-shot injections for essentially complete 
combustion, cutting emission of particulate matter and NOx. In other applications the present 
invention may permit control of ultra-short opening and closing of the primary solenoid and short 
controllable dwell interval between two impulses (or a series of impulses). In other words, under 
the present invention the dynamic time series may become very close to electromagnetic wave 
forms indicated by an electric signal output from the actuators. 

Referring now to Fig. 9 (with its x-axis coordinate setup), it is seen that at the start of 
injection 0 < / < r , while the needle moves upwards, a force accelerating the needle valve 
with mass m is superposed by: magnetic force FMag induced by an energized solenoid 
(primary coil), elastic force Fei produced by a compressed spring, gravity force Fgr due to 
universal Earth gravitation (9.98 m/s^) and side friction force Ffr because of contact of the 
needle surface to a thin fuel layer occurring in a high pressure fuel passage: 
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F^^ = 5//sin(0«) = //,MoM/ = u,u,PN 



Fgr=mg 



(2) 
(3) 
(4) 
(5) 



where B is magnetic flux density (induction), is relative permeability of ferromagnetic 
iron, Uq= 1.257*10'^ H/m is magnetic field constant, / is coil (solenoid) length, / is 

current supplied to coil, N is number of turns on coil, k is spring constant according to 
Hooke's law, is initial spring compression, and q^^^ is friction coefficient under laminar 

conditions (turbulent component of the friction force is neglected due to very thin layer in 
the fuel passage resulting in low Re-number). 



Temporal transition conditions are: 



r = 0: 



/ = Q{A\x = ^^[m\x' = Q\mls] 



(6) 



t = r: 



I^I^[Alx^{^,+^)[m■\ 



(7) 



In general, an exponential law presents transient time dependent current: 



1=1 Jit) 



(8) 



Now, eq. (1) can be rewritten in the form of: 



m 



X X- 

m m 



m 



or 



(9) 



x' + af,x' + a^fX = a„^JJ (t)-a 



sys 

n-linear linear 



The above implies that a solution of this second-order non-homogeneous ordinary 
differential eq. (9) will be obtained using superposition of two exponent type ftinctions 
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x(t) = x^ (t) + X2 (0 of the arguments dependent on time t and amplification factors y , so they have a 
transient oscillatory nature during start-up of transition, with respect to linear and non-linear parts 
on the right hand. The first function regarding linear part of the eq. (9) has a generalized form as 

x,(t)^Aoe^' (10) 

Using derivatives of x' and x'^from the function in the eq. (9), the linear part 
becomes in form of: 

Ao(A'+«/.A+«e/)^'^' (11) 

At very beginning, vy^hen / -> 0, this expression is transferred to a quadratic equilibrium: 

fi^ + a^^fi, + (a,, + a^, /A,) = 0 (12) 

which can be resolved with respect to variable JJ^ , i.e. basic frequency of oscillation: 



In general, there are three classes of solution depending on the sign of square root in eq. 
(13). However, in the case of solenoids applied to move a needle inside of a high pressure 
fuel barrel, for example, the friction force is negligibly small versus elastic and gravity forces 
^fr « 4(«e/ + ^sys) » ^^c solutiott to basic frequency /? can be rewritten as: 



A = ±^a^i+a,^JA, = ±10), (14) 

and general solution Xi(t) for the upward lifting dynamics at the start of injection is: 
X, (0 = Aoe*"^' = Ao[cos(fi>iO ± is\n(a),t)] (15) 

The second function regarding non-linear part of the eq. (9) has the same generalized 
form as: 

x,(t) = r,e^^' (16) 

Taking derivatives of x' and x" from the fUnction X2(t) in the eq. (9), one can obtain 
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equilibrium of: 



(fil +a,J, +a„)r,e^' =a^llfit) (17) 

Given an electric circuit of solenoid composed of an inductor with inductance L and a 
resistor with resistance R in series connection, the Kirchhoff loop rule requires that the 
sum of the changes in potential around the circuit must be zero, so: 

Z— +/i? = 0 (18) 
The solution for this eq. (18) is: 

(19) 

The magnetic field of a current-carrying conductor or a coil changes with the conductor 
current. A voltage proportional to the change in current is induced in the conductor itself 
and counteracts the current change producing it. Therefore, for the self-induction, eq. (18) 
is transformed to: 

-i— + /i? = 0 (19.1) 
dt 

which solution is: 

I = lJ^ (20) 

Now, assuming only one solenoid or coil forcing the needle upward, which current is 
described by eq. (19), one can rewrite (17) as: 

{p\ +a,,)Yy^' = a^^lli (21) 

from which the solution can be found using equality of constant and time dependent 
parts: 

{Pl^a^rP2^0Le,)y2=CCn.A 0-2) 
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e^^'=e ' (23) 

and general solution, expressed by eq. (16), assuming negligibility of friction force versus 
magnetic and elastic forces, becomes: 



^2(0 = 



-os^^ (24) 



— + « 



where sign reflects start up (switch-on) of the solenoid and "-" reflects switch off 
of the solenoid, 0^21 a transient frequency determined time response, k is amplification 
factor due to combination of the injector and solenoid construction parameters, and is 

a current level which is limited because resistance heat-cooling balance suffering bum 
damage. This second lift component (t) is much greater than (/) while the solenoid 

of injector (or of an actuator) is energized. The time response is limited by all three factors 
indicated in eq. (24) and for a given injector/solenoid configuration can be controllable only 
through possible control (increase) of transient frequency a>2i . 

Now, assuming that at the transient moments the current applied to a primary coil 
characterized by , and 6)2^ , is generated by a remote (not installed physically on 
the same injector or actuator housing) solenoid characterized by ^2 » hi ^22 > which is also 
just energized or de-energized (opened or closed). Transmission of the self-induction 
transient current from secondary solenoid to the first coil will generate a very special sharply 
shaped current that can be performed by super-exponential "I-Function": 



(25) 



This function operates as a modulation function /(/) in eq. (17), i.e., it implies a speed of 

dynamic influencing directly on transient frequency (or time response) of the primary "physically" 

installed solenoid. Some basic features of the I-Function and its first order derivative are shown in 

Fig. 10. As seen in this Fig., the maximum peak phase of the current is gradually shifted upon a 

R / 

magnitude of 6^22 (in other words by a factor VC of the secondary coil) while the peak amplitude 

/ ^2 

is dependent on (o^y (in other words by a factor ). The transition period is also controllable 

depending upon the ratio between co^x and 6)22 • The higher magnitude of this ratio determines the 
shorter transition. 

The same ratio factor controls the speed of lifting indicated by the first order derivatives: 
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the higher ratio y reflects more rapid speeding of the needle lift. The turnover points in the 

/ ^22 

bottom plot of Fig. 10 indicate that rapid "one^peak" acceleration is achieved at higher ratio values. 
The lower ratio may reflect a series of acceleration peaks. Of note, the secondary solenoid may be 
presented by a nonphysically installed remote coil. It can be also coded as a signal (e.g., a digital 
signal) and, using a D/A converter, for example, supplied to the primary coil. An illustrative 

secondary and primary coil configuration may utilize a highest ratio of that excludes a 

longer transition and makes possible to induce strong magnetic flux in the primary coil within 
shortest time permitting a long time of heat dissipation (e.g., the shortest transient induction duty 
cycle permitting afterwards to run ultra-shot multi-injection cycle per each injection stroke). 
Criteria to select operation parameters of the coils are determined by the momentum 
equations: 

{Pi ^a^^P, ^a,)y,e^^^ =a^^jll,{t) (26) 
which implies that: 

{Pl^a^J,+a,,)Y,^a^jl (27) 
p^N^dl^ (28) 

The first equation (27) determines construction of the primary coil in terms of inductance 
i, and time response VC . The second equation (28), the rapid speeding, permits to 

calculate ratio of which is used for deduction of the secondary coil properties: inductance 

/ ^22 

I2 and time response 2/ or take out the input signals to a secondary solenoid digital (electronic) 
model. 

Referring now to Fig. 1 1 A, one example of an electric circuit incorporated secondary coil 
(which example is intended to be illustrative and not restrictive) is shown. More particularly, Fig. 
1 1 A shows a simple inductive pre- and post secondary inductor circuit (e.g., for a fuel injection 
system) and Fig. 11 B shows two associated timing scenarios. In these Figs, the secondary inductor 
or secondary coil (SC) is designed to create a fuel injector driver, which uses one or two secondary 
inductors to enhance injector performance. Of note, this equipment may generate much higher 
voltages than normal fiiel injector drivers, which may break the injector's dialectic insulation and/or 
can cause injury to the unwary operator. Therefore, critical parameters may first be simulated using 
code (e.g., the code described below). In addition, although quicker fiiel injector currents are 
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expected, there is no guarantee as to the physical speed or change in speed of the injector. 
Therefore, each new model may be verified using specially developed test equipment. Later on, one 
can find a description of the test procedures with regard to the fiiel injectors for internal combustion 
engines 

In any case, the circuit in Fig. 1 1 A may operate as follows: 

Before the injector solenoid with inductance LI is fired the secondary inductors, L2 and L3 
will be pre-charged. Both transistors Tl and T2 are turned on at this time. 

Transistor Tl is turned off when injection is desired. 

The current, pre-charged on the secondary L2, generates a high voltage that drives the 
injector inductor, i.e. primary coil ("PM"). 

Afterwards, the current stabilizes to maintain the valve open. 

Turning off transistor T2 leaves currents in the injector (LI) and inductor (L3) competes 
causing much higher voltages at TP2. The competing currents will also terminate the 
injector current quicker. 

Of note, the circuit schematic of Fig. 1 1 A represents system basics generically, not 
specifically to the final circuit related to specific injector and/or other type of actuators. For 
instance, the secondary inductors may be varied and additional resistance may be added for steady 
state operations. The main driving transistors may also require their own drivers. The charge time is 
easily controlled through the charge time of L2. The Rl is the resistance added in the driver. That 
resistance is essentially only to safeguard the circuit. If the L2 charges too long the circuit may bum 
up. In final configuration, the ECU of the vehicle may protect the final circuit. The transistors are 
treated as switches, so they are ignored for the purposes of the simulation code discussed below. 
Since Tl is off and T2 is on, for the simulation program it is necessary to consider the current 
stream going from parallel CI -LI loop farther through chain of injector components R3-L1-R4 to 
transistor T2. The T3 is in the event that a ftinction generator could not drive the Tl transistor. The 
Tl transistor only has an amplification of about 12, hence it takes almost 1 amp for the transistor to 
drive 10 amps. To get supercharging of the secondary coil the electric circuit may need to be 
changed in such way that the secondary coil is connected to the primary injector coil skipping over 
the control resistor (in Fig.l lA the connection of L2 is going directly to LI skipping over R3). 
One may need to drive transistors Tl and T2 through Rl and R2 respectively with a control 
device capable of 1 amp power supply. The values are dependent on voltages. Care may need to be 
applied in selecting the proper transistors (although MOSFETs are typically cheaper and easier to 
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design with, practical experience shows that a good Bipolar may survive test more reliably). 
Accordingly, while various circuit parameters may be changed as desired and/or dictated by 
application, it is understood that such changes are readily within the reach of those of ordinary skill 
in the art in light of the present disclosure. 

Referring now to code for the calculation of the secondary coil charging time (an example 
of which code is described below), it is noted that such code may compute a minimum time needed 
to charge a secondary coil for generating an I-Function like shaped current depending on 
inductance and resistance characteristics of the primary and secondary coils as well as initial 
current and voltage values applied to the capacitor and the coils. Direction of the current through 
secondary coil L2i and Llj as well as voltage onto the capacitor Cv are schematically indicated in 
Fig. 1 1 A. The calculation is based on basic current and voltage equations applied to a capacitor and 
an inductor: 

dV 

Ic=C^ (29) 

K=C^ (30) 
at 

where Fand / are time dependent variables. The change in voltage on the capacitor is: 

dc,=^^^lzHLdt (31) 
C 

In addition, the voltages associated with resistances of secondary R2 and primary 
Rl coils are: 

Rly^Ll^Rl (32) 
R\y=L\,R\ (33) 

From Fig. 1 1 A one can write the voltages balance on secondary L2v and primary LI v coils 

as: 

L2y=V^,,,^-Rly^C, (34) 

L\y^Cy-R\y (35) 



Therefore, according to equations (29) and (30), the changes in current through secondary 



and primary coils can be derived to: 




(36) 




(37) 



Turning now to a specific example of computer code for determining various parameters 
associated with the present invention (which example is intended to be illustrative and not 
restrictive), the following code may be used: 

program secondary solenoid 

c +12V — L2 - R2 -o- LI Rl -- Gnd 

c +.+ I 

c C 

c I 
c Gnd 



c Ic = Cdv/dt->dv = Ic/C*dt 
c Vi = Ldi/dt-o>di = Vi/L*dt 

real L2,L1,R2, Rl 

real L2i, Lli, L2v, Llv, R2v, Rlv 

real t, dt 

real C, Cv, Vin 

integer i 

c 

c input basic parameters 

open (4,file='Input_Electric.dat') 

read (4,'(a80)')dummy 

read (4,*) L2 

read (4,'(a80)')dummy 

read (4,*) R2 

read (4,'(a80)')dummy 

read (4,*) LI 

read (4,'(a80)')dummy 

read (4,*) Rl 

read (4,'(a80)')dummy 



c 
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read (4,*) C 

read (4,'(a80)')dummy 

read (4,*) Vin 

read (4,'(a80)')dummy 

read (4,*) L2i 

read (4,'(a80)')dummy 

read (4,*) LI i 

read (4,'(a80)')dummy 

read (4,*) R2v 

read (4,'(a80)')dummy 

read(4,*)Rlv 

read (4,'(a80)')dummy 

read (4,*) Cv 

read (4,'(a80)')dummy 

read (4,*) t 

read (4,'(a80)')dummy 
read (4,*) dt 
read (4,'(a80)')dummy 
read (4,*) Nt 
close (4) 



open (10,file='AI IData.dat') 
write (10,*) 'L2', L2*le3, ' [mH]' 
write (10,*) 'R2', R2, ' [Ohm]' 
write (10,*) "LI", Ll*le3, ' [mH]' 
write (10,*) 'Rl'.Rl,' [Ohm]' 
write (10,*) 'C=', C*le6, ' [uF]' 
write (10,*) 'Vin=', Vin, ' [V]* 
write (10,*) •L2i', L2i, ' [A]' 
write (10,*) 'R2v', R2v, ' [V]' 
write (10,*) 'Lli', Lli, ' [A]' 
write (10,*) 'Rlv', Rlv, ' [V]' 
write (10,*) 'Output Data:' 

write (10,*) 'L2 charge time=', L2i*L2A^in/le-6, ' [ us]' 
write (10,*) •t[us] Cv[V] L2i[A] Lli[A]' 
do i=l,Nt 

Cv = Cv + (L2i-LIi)/C*dt 
if(Cv.le.-1.4)Cv=-1.4 
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R2v = L2i * R2 
Rlv = Lli*Rl 

L2v = Vin - R2v - Cv 
Llv = Cv-Rlv 
L2i = L2i + L2v / L2 * dt 
Lli = Lli + Llv/Ll * dt 
write (10,89) t*le6, Cv, L2i, Lli 
89 format (f5.1, 2x, f6.1, 2x, f5.1, 2x, fS.l) 
t = t + dt 
enddo 
close(lO) 
stop 

end 

Input Data File 

L2 is inductance of secondary solenoid, [H] 
0.000209 

R2 is resistance of secondary solenoid, [Ohm] 

0.5 

LI is inductance of primary (injector) solenoid, [H] 
0.0005 

Rl is resistance of secondary solenoid, [Ohm] 
20.0 

C is capacity, [F] 
0.33e-6 

Vin is supply voltage, [V] 
24.0 

L2i is initial current through secondary solenoid, [A] 
8.0 

LI is initial current through primary (injector) solenoid, [H] 
0.0 

R2v is initial votage applied on secondary solenoid, [V] 
0.0 

Rlv is initial votage applied on primary (injector) solenoid, [V] 
0.0 

Cv is Initial volage on capacitor, [V] 
0.0 

t Is initial time, [s] 
0.0 
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dt is time increment, [s] 
2.0e-7 

Nt is number for timing, [-] 

1200 

M is number for data print control 
10 

Output Data File 

L2 0.209000006 [mH] 
R2 0.500000000 [Ohm] 
LI 5.00000000 [mH] 
Rl 1.29999995 [Ohm] 
C= 0.330000013 [uF] 
Vin= 24.0000000 [V] 
L2i 8.00000000 [A] 
R2v O.OOOOOOOOE+00 [V] 
LliO.OOOOOOOOE+00[A] 
RlvO.00OOO0O0E+OO[V] 
Output Data: 

L2 charge time= 69.6666718 [ us] 

t[us] Cv[V] L2i[A]Lli[A] 

0.0 0.0 8.0 0.0 

2.0 53.3 7.9 0.0 

4.0 99.8 7.3 0.0 

6.0 141.4 6.4 0.1 

8.0 175.7 5.0 0.2 

10.0 200.6 3.4 0.2 

12.0 214.8 1.6 0.3 

14.0 217.5 -0.2 0.4 

16.0 208.4-2.0 0.5 

18.0 188.2 -3.7 0.6 

20.0 158.3 -5.1 0.6 

22.0 120.3 -6.2 0.7 

24.0 76.6 -6.8 0.7 

26.0 30.0-7.0 0.8 

28.0-1.4 -6.9 0.8 

30.0-1.4-6.6 0.8 

32.0-1.4-6.3 0.8 

34.0-1.4-6.0 0.8 



36.0-1.4 -5.8 0.7 
38.0-1.4 -5.5 0.7 
40.0-1.4 -5.2 0.7 
42.0-1.4 -5.0 0.7 
44.0-1.4 -4.7 0.7 
46.0-1.4-4.4 0.7 
48.0-1.4-4.2 0.7 
50.0-1.4 -3.9 0.7 
52.0-1.4 -3.6 0.7 
54.0-1.4 -3.4 0.7 
56.0-1.4 -3.1 0.7 
58.0-1.4 -2.9 0.7 
60.0-1.4-2.6 0.7 
62.0-1.4-2.4 0.7 
64.0-1.4-2.1 0.7 
66.0-1.4-1.8 0.7 
68.0-1.4-1.6 0.7 
70.0-1.4-1.3 0.7 
72.0-1.4-1.1 0.7 
74.0-1.4 -0.8 0.7 
76.0-1.4 -0.6 0.7 
78.0-1.4 -0.4 0.7 
80.0-1.4 -0.1 0.7 
82.0-1.4 0.1 0.7 
84.0-1.4 0.4 0.7 
86.0-1.4 0.6 0.7 
88.0-1.2 0.9 0.7 
90.0 0.2 1.1 0.7 
92.0 3.0 1.3 0.7 
94.0 6.9 1.5 0.7 
96.0 11.8 1.6 0.7 
98.0 17.3 1.7 0.7 
100.0 23.1 1.7 0.7 
102.0 28.9 1.7 0.8 
104.0 34.2 1.6 0.8 
106.0 38.9 1.5 0.8 
108.042.5 1.3 0.8 
110.0 45.0 1.1 0.8 



112.0 46.1 0.9 0.8 
114.0 45.7 0.7 0.8 
116.0 44.0 0.5 0.9 
118.041.1 0.3 0.9 
120.0 37.0 0.1 0.9 
122.0 32.1 0.0 0.9 
124.0 26.7 0.0 0.9 
126.0 21.0 0.0 0.9 
128.0 15.5 0.1 0.9 
130.0 10.4 0.2 0.9 
132.0 6.1 0.3 0.9 
134.0 2.8 0.5 0.9 
136.0 0.7 0.7 0.9 
138.0-0.1 0.9 0.9 
140.0 0.5 1.2 0.9 
142.0 2.4 1.4 0.9 
144.0 5.5 1.6 0.9 
146.0 9.7 1.7 0.9 
148.0 14.6 1.8 1.0 
150.0 19.9 1.9 1.0 
152.0 25.4 1.9 1.0 
154.0 30.7 1.8 1.0 
156.0 35.5 1.7 1.0 
158.0 39.5 1.6 1.0 
160.0 42.5 1.4 1.0 
162.0 44.2 1.2 1.0 
164.0 44.7 1.0 1.1 
166.0 43.8 0.8 1.1 
168.0 41.6 0.6 1.1 
170.0 38.4 0.5 1.1 
172.0 34.1 0.4 1.1 
174.0 29.3 0.3 1.1 
176.0 24.0 0.3 1.1 
178.0 18.7 0.3 1.1 
180.0 13.6 0.4 1.2 
182.0 9.1 0.5 1.2 
184.0 5.4 0.6 1.2 
186.0 2.8 0.8 1.2 



188.0 1.3 I.O 1.2 

190.0 1.2 1.3 1.2 
192.0 2.3 1.5 1.2 
194.0 4.6 1.6 1.2 
196.0 8.0 1.8 1.2 
198.0 12.3 1.9 1.2 
200.0 17.1 2.0 1.2 
202.0 22.3 2.0 1.2 
204.0 27.4 2.0 1.2 
206.0 32.2 1.9 1.2 
208.0 36.4 1.8 1.2 
210.0 39.8 1.7 1.2 
212.0 42.1 1.5 1.2 
214.0 43.2 1.3 1.3 
216.0 43.1 1.1 1.3 
218.041.7 1.0 1.3 
220.0 39.2 0.8 1.3 
222.0 35.7 0.7 1.3 
224.0 31.4 0.6 1.3 
226.0 26.6 0.5 1.3 
228.0 21.6 0.5 1.4 
230.0 16.7 0.6 1.4 
232.0 12.1 0.7 1.4 
234.0 8.1 0.8 1.4 
236.0 5.1 1.0 1.4 
238.0 3.1 1.1 1.4 



Referring now to secondary coil charging scenarios and electric wave forms, it is noted that 
at least two different charge-timing scenarios may be applied. In one, the secondary coil SC 
is charged (e.g., from zero to a few thousands of microseconds) essentially simultaneously with the 
injection duration signal applied to the primary coil (PC), in other words, essentially 
simultaneously with the primary coil. As seen in the bottom part of Fig. 1 IB, the charging period of 
the SC is controlled by the transistor Tl and indicated by triggering impulse AB. Closing, opening, 
and closing of the PC is controlled through transistor T2. Impulse CD at the transistor indicates 
injection duration pulse. This scenario is called "simultaneous charge". 

In the second scenario, the SC is charged first and afterwards a signal is applied to the PC. 
In Fig. 1 IB this is shown as series of triggering impulses AB at Tl and CD at T2. This scenario is 
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called "pre-charge" (there is another scenario when the SC starts charging and during this phase, 
after some delay, the PC also starts its duty cycle (injection duration signal at T2); this mixed 
charging scenario is called "shifted charge). 

Fig. 12A illustrates typical waveforms for simultaneous charge of the SC and Fig. 12B 
illustrates typical waveforms for the pre-charge of the SC. Because of the inductance of the SC in 
the circuit and connection of L2 in series with LI, in both cases the charging of the PC starts with 
delay essentially equal to the time at which the SC is charged. However, the waveforms obtained 
from a tested injector are different. 

Under simultaneous charge, the diagram in Fig. 12A, the magnetic energy accumulated into 
the SC transfers rapidly and at higher level of amplitude. Two phase-separated spikes are observed. 
The first spike shows start of the SC charge. The second spike indicates startup of the PC operation 
(injection duration). This regime is very important for injection and combustion control (e.g., over 
diesel engines). It permits the split of the whole injection cycle per each stroke in multi-shot ultra- 
shot injection series (e.g., pilot injection and series main injection). This allows, as seen in Fig. 13, 
the transfer of a diesel stratified diffusion flame structure into a "Christmas-like" structure with 
multi access of air into the diffusion flame boundaries (resulting in more complete combustion at 
any given rate of the fuel; increase in fiiel economy; and a cut-off emission of particulate matter 
and NOx). 

Referring once again to Fig. 12B, it is seen that this diagram relates to the "pre-charge" 
case. The first spike indicates charging of the SC and in "cascade" the second spike shows charging 
of the PC and startup the injection. At the transition moment one can see a small "zigzag" type 
oscillation which indicates that the PC is rapidly interfered with magnetic flux of the SC. This 
regime is particularly applicable for gasoline engines (especially for direct injection gasoline 
engines where the spray structure is stratified). Rapid opening of the valve permits the spray to 
reach fine quality within very short time fi*action. If the injector has swirl nozzle exit, this technique 
permits control of swirl speed (rotational speed) that results in a fine spray essentially immediately 
after fijel jet breaks up into the spray. The same case is important for the diesel engines at the 
moment when one needs to organize multi-shot injection, described above (e.g., a main injection 
with well controlled dwell intervals between injection shots). 

Referring now to verification of injection system operation (e.g., speed), it is noted (as 
mentioned before) that there is no guarantee regarding the timing response of the whole injector 
system (i.e., even if the electric output signal from the fiiel injector coupled with the SC controller 
indicates fast response). Direct applications of a secondary solenoid (SC) in automotive field are 
typically related to diesel and direct injection gasoline engines where a stratified charge of fuel 
mixed with tumbled or vortex airflow determines the quality of combustion and its completeness. 
The spraying of fuel typically ends immediately after dropping down the pressure in the 
accumulation injector chamber (or high-pressure gallery). In other words, the closing timing on the 
valve is a quite rapid process because propagation of the pressure waves with sound speed brakes 
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the spray even before the mechanical sealing of needle at the nozzle exit occurs. So, in one 
embodiment the concentration is on the valve opening process. 

In this regard, the focus may be placed on injection shot duration ("ISD") with controllable 
rise time and holding time and the dwell interval ("DI") between the shots. In one example (which 
example is intended to be illustrative and not restrictive) relating to common rail diesel injectors 
(e.g., a Bosch system) the ISD is matched at a few tens of microseconds (comparable with "fuel jet 
break-up time) and DI is matched at a few hundreds of microseconds (limiting to oxidation cycle 
per single shot to keep diffusion flame around the core spray). 

The pilot injection and main injection may to be split into a multi-shot injection series. In 
DI gasoline engines these requirements may be different; instead, it may be necessary to have only 
one -100 ms shot phased properly to the igniting moment. To make a robust and simple 
verification of SC impact and operation, one may have an injection system with initially 
controllable injection period (T) and injection duration (tau). 

A configuration of a system for managing the injection flow according to an embodiment of 
the present invention is shown in Fig. 14. A control signal from a sensor (or any available feedback 
line) is fed to the ECU receiving the signal from all sensors on the engine board and transmitting 
control signals to the execution parts of the engine. The ECU output also manages the injector 
primary coil (PC) in terms of current and/or voltage applied onto the PC and depending on the 
engine run regime produces a current and/or voltage applied onto the secondary coil SC. The SC 
generates an I-Function like current and the injector rapidly starts to operate (rapid opening of the 
valve due to magnetic flux). 

In order to help ensure that rapid opening of the valve actually takes place (not only electric 
wave front obviously seen on oscilloscope), the control measurement may be done using the LDV 
Instantaneous Flow Rate Measurement Stand described in applicant's pending U.S. Patent 
Application 20020014224, published February 7, 2002. 

For a demonstration of this rapid response even at low injection pressure, the inventor has 
built up a test cell, which simulates the injection system depicted in Fig. 14 and described above. 
The test cell is depicted in Fig. 15 and composes four sub-systems: 

• The injection system is represented by a fuel tank pressurized by inert nitrogen gas. The fuel 
delivering line is connected to a measurement intersection in which a capillary quartz pipe is 
installed. The measurement intersection is constructed to operate at both steady state and 
oscillatory fuel flow under high injection pressures generated in diesel injection systems. The 
metal intersection itself is mounted in heavy metal frame with 3D alignment and adjustment 
mechanics. The outlet of the measurement intersection is flexible to mount essentially any type 
of fuel injector. 

• A Laser Doppler Anemometer ("LDA") of Dantec/Invent Measurement Technology GmbH is 
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used to measure centerline velocity into the fuel flow oscillating in the quartz pipe. The LDA 
consists of the Transmitting and Photo-receiving Optics, the Ion Laser coupled to the Fiber 
Transmission units, the Fiber PDA 58N70 Detector Units, the Multi PDA 58N80 Signal 
Processor and the Dantec 3D Traverse. An LDA signal can be observed using the Hewlett 
Packard Infmium 500 MHz IGsa/s Oscilloscope. To monitor cyclically operating injection 
flow, the Cyclic Phenomena Dantec software is used to process and treat the output results. 
Angular encoding signal is provided from a Waveform Generator (e.g., the same one which 
controls injection duty cycle). The system measures forward and reversed velocity due to the 
Bragg Cells in the transmitting optics. The main parameters used for the demonstration 
measurements are: 

o Optical probe 77x77x945 mm 

o Fringe spacing 3.15 mm 

o Frequency shift 40 MHz 

o Cyclic length 360 degree 

o Phase averaging bins 360 

• The injector driver system starts from the Agilent 33 120 A 1 5 MHz Function/Arbitrary 
Waveform Generator which precisely controls TTL signal frequency. The Stanford Research 
System, Inc. Model DO 535 Four Channel Digital Delay / Pulse Generator has 8 input/output 
ports that used to adjust various delays with respect to initially generated TTL trigger impulse 
waveform. Particularly, AB and CD ports are used to control charging time of secondary coil 
by transistor Tl and injection duration of injector primary coil by transistor T2, respectively, A 
Regular automotive battery of 12 V is used as the DC power supply. The output voltage firom 
the secondary coil driver is directly connected to the test injector. The injector plug unit has 
input/output ports, so the output signal is observed at the Tektronix 2221 100 MHz Digital 
Storage Oscilloscope. 

• To verify accuracy of the LDA flow rate measurements, the injected mass time series are 
recorded using the A&D Company, Ltd. GX-4000 Multi-Functional Balance (simultaneously 
with the LDA time series). Measurements in steady state and oscillatory flows shows that in 
laminar flow accuracy lays within 1.1 %, in turbulent flow it comes within 2.3 %. 

In the above example all demonstration measurements were conducted under pressure of 
7.3 atm (105.85 psi) at the injection frequency of 50 Hz (20 ms cycle period). Two different 
charge-timing scenarios were applied. Firstly, SC coil was charged from zero to 2000 microseconds 
and afterwards the primary solenoid coil (PC) was opened. Injection duration in this particular 
example was the same for all measurements of 15 ms. Secondly, the secondary coil was charged 
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from zero to 2000 microseconds simultaneously with the injection duration signal applied to the 
primary coil. Injection duration was setup at 3 and 5 ms, at each case a number of the instantaneous 
flow rate time series were measured. 

Referring now to computer code for operating on each centerline velocity time series 
associated with the present invention, one example of such computer code (which example is 
intended to be illustrative and not restrictive) may be as follows (of note, this program reconstructs 
the measurement data into instantaneous series of volumetric/mass flow rate, pressure gradient and 
integrated (or accumulated) fuel mass within each injection cycle): 

c For Turbulent Flows 

program FlowRate_MSU_07 
external bessjO,bessj 1 
complex bessjO, bessjl 
complex i 

real tint, M_mean, M beg, M_per, M Jnt 
character*2 Al, fiiame*12 
complex Q(4096), C(4096), P(4096) 

real U(8192), UB(8192), U_t(8192), ph(8192), U_cor(150,150) 

real Qcor(8192), P^Z(8192), Q_u(8192), Massjnt(8192) 

integer Nexp, 1, j, NP, NR 

real nue, rho, TO, R, tau, k, d_tph 

c 

c input basic parameters 

open (4,file='Input_Fuel_BKM.datO 

read (4,'(a80)*)dummy 

read (4,*) TO 

read (4,'(a80)0dummy 

read (4,*) nue 

read (4,'(a80)')dummy 

read (4,*) rho 

read (4;(a80)')dummy 

read (4,*) R 

read (4,'(a80)')dummy 

read (4,*) tau 

read (4,'(a80)Odummy 

read (4,*) k 

read (4,'(a80)')dummy 
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read(4,*)NR 
read (4;(a80)')dummy 
read(4,*)NP 
close (4) 



fO=l./TO 

i = (0.,l.) 

pi = 4.*atan(l.) 

w0 = 2*pi*f0 

TeO = R*sqrt(wO/nue) 

c 

c input array of the measured velocity series 

c within the period using "Ivr" software, TO is equal 720 degree 

open (5,file='ldv.dat') 

1=0 

10 1=1+1 

read(5,*,end=12) nn, ph(l), nl, u(l), rms 
c REVERSED Measurement! 
u(l)=(-iru(l) 
goto 10 
12 continue 
close(5) 

write (*,*) 'experimental data file have been read' 
Tint= TO 
Nexp= 1-1 



c avarage parameters obtained from direct velocity 
c time-series measurement 

doof-0. 

do 1 = l,Nexp 

doof = doof + u(l) 
Q_u(l)= u(l)*pi*R*R/2. 

enddo 
c mean of velocity 

U_mean = doof/float(Nexp) 
c mean of mass rate 

M^beg = U_mean*pi*R*R*0.697*rho 
c mean of mass per one statistical cycle 
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Mjer = M_beg*Tint/1000 

c 

c Fourier transform and its inverse 

c with respect to equidistant time-phases ph(l) 

call ffi (u,C,Nexp) 

call ffs (ub,C,Nexp) 

open (6,file='check.dat') 

do j= l,Nexp 

write (6,*) phO),uOXubG) 

enddo 

close (6) 

write (*,*) 'passed Fourier transform and its inverse' 

Q=:^==== = === =: 7 -^ r = 

c complex components of pressure gradient 
c normalized by density rho 

open(66, file-prgr_comp.dat') 

P(l)= C(l)*2.*nue/(R*R) 

write(66,*) real(P(l)), imag(P(l)) 

doj= 2,Nexp/2+l 

Ten = R*sqrt(G-l)*wO/nue) 

PG)= CC)*0-l)*wO*i/(l.-l./bessjO(i**L5*Ten)) 

write (66,*) reai(P(j)), imag(PG)) 

enddo 

write (*,*) 'normal.compl.component of press.gradient' 
c================================— ===== 

c computing the theoretical velocity time-series 
c on a pipe axis 

open (7,file='theory.dat') 

do ln= 1, 100 

Uj(ln)- P(l)*R*Ry(4.*nue) 

tph= float(ln)/float(Nexp)*2,*pi 

doj= 2,Nexp/2+l 

Ten = R*sqrt(0-I)*wO/nue) 
wn= wO*(j-l) 

U_t(ln)= Real( U_t(ln)+ PO)*i*cexp(i*tph*G-l))/wn* 
& (1 ./(bessjOG* * 1 .5*Ten))-l .)) 

enddo 

write (7,*) ph(ln), ub(ln), U_t(ln) 
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write (*,*) ph(ln), ub(ln), Uj(ln) 

enddo 

close (7) 



c complex component of flow rate 
c open (77,file='compl_FR.dat') 

Q(l)= 0.697*P(l)*pi*R**4/(4 *nue) 
c write (77,*) Q(l) 

doj=2,Nexp/2H-l 

Ten = R*sqrt((j-l)*wO/nue) 

Q(j)= 0.697*PG)*pi*R*R*i/(wO*G-l))* 
& (4.*i**0.5*bessjl(i**1.5*Ten)/(Ten*bessjO(i**1.5*Ten))-2.) 
c exponensial oscillation is given below 

write (*,*) QO) 

enddo 



c computing of flow rate time-series 
c and avarage parameters 
Q_int= 0. 

d_tph= TO/float(Nexp) 
do ln= 1 ,Nexp 
Qcor(ln)= Q(l) 

tph= float(ln)/float(Nexp)*2.*pi 
doj= 2,Nexp/2+l 

Qcor(ln)=real(Qcor(ln)+Q(j)*cexp(i*tph*G-l))) 
enddo 

QJnt= CLint+Qcor(ln) 
Mass_int(ln)= Q_int*rho*d_tph 

enddo 

c mean of mass per one period 
Mint = Q_int/float(Nexp)*rho 
M^mean = Real(Q(l))*rho 

write (*,*) 'flow rate was integrated' 



c computing of pressure gradient 
do ln=l,Nexp 
P_Z(ln)= P(l) 

tph= float(ln)/float(Nexp)*2.*pi 
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doj= 2,Nexp/2+l 
P_^Z(ln)= P_Z(ln) + Pa)*cexp(i*tph*0-1)) 

enddo 
P_Z(ln)=-rho*P_Z(ln) 

enddo 

write (*,*) 'pressure gradient was computed' 



open (10,file=' AllData.dat') 

write (10,*) •CA[deg] U[m/s] Vj[ml/s] P_z[MPa/m] Massjnt[g] ' 
do ln= l,Nexp 

write (10,89) ph(ln), u(ln), Qcor(ln)*1.0e6, P_z(ln)/K0e6, 
&MassJnt(ln) 

89 format (f6.1, 2x, f7.3, 2x, f73, 2x, f9.5, 2x, f8.5) 
enddo 
close(lO) 

open (1 1, file-result.dat') 

write(l l,*)'Injection cycle TO:^TO,'[ms]' 

write(l l,*)'Mean velocity U mean :',U_mean,'[m/s]' 

write(l 1,*)'MR: di vel int M^beg :',M_beg,'[kg/s]' 

write(l l,*yM/cycle: si vel int M_per:',Mj)er,'[kg]' 

write(l l,*)'Integrated mass flowrate M_int:',MJnt,'[kg/s]' 

write(l 1,*)' *Mass: the first Fourier term:*,M_mean,'kg/s]' 

close(ll) 

stop 
end 

c= 

complex function bessjO(x) 

external summe 
complex X 

complex summe,bess 
integer j 
bess=(l„0.) 

dorl,12 

bess=bess + summe(xj) 

enddo 

bessjO=bess 

return 
end 
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c 

complex function summe(z,n) 

integer n 

real prod 

complex z 

5 

prod=L 
doj= l,n 

prod= prod*float(j) 
enddo 

prod= prod*prod*((-l)**n) 

summe= (0.25*z*z)**float(n)/cmplx(prod) 

return 

end 

c 

complex function bessjl(x) 
external summel 
complex X 

complex summel, bess 
bess=(0.,0.) 

doj=l,12 

bess= bess +summel(xj) 

enddo 

bessjl=bess 

return 

end 

c 

complex function summel (z,n) 

integer n 

real prod 

complex z 

prod=l. 

do J=I,n 

prod=prod*float(j) 
enddo 

prod =((-0.25)**n)*2.*float(n)/(prod*prod) 
summel =prod*(z* *float(2*n- 1 )) 
return 
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end 

C======== — = •" —======== 

subroutine ffl(X,C,N) 
integer N 

complex C(4098), pin 
real X(8192) 
do i=0,N/2 

pin - (O.,10*(8.*atan(L)*dble(i)/dble(N)) 

C(i+1H0.,0.) 

6 

dorUN 

C(i+l)=C(i+l)+dcmplx(X(j))*CDEXP(pin*dcmplx(-j)) 
enddo 

C(i+l)=C(i+l)*dcmplx(2./dble(N)) 

enddo 

return 

end 

C ===== = ===== == = 1 = =.-■ :=====:==== 

subroutine fFs(X,C,N) 
integer N 

complex C(4098), argum 
real x(8I92) 
do i=l,N 

argum = (0.,l.)*(8.dO*atan(l.)*dble(i)/dble(N)) 

x(i) = dble(C(l)*0.5) 

doj=l,N/2 

x(i) = x(i) + dble(C(j+l)*cexp(argum*j)) 

enddo 

enddo 

return 

end 

Three different SC charging techniques are depicted in Figs. 16A and 16B. All the data in 
these Figs. 16A and 16B were measured under the same conditions: injection frequency 50 Hz, 
injection pressure 7.3 atm and SC charging time 2.0 ms. Fig. 16A shows instantaneous volumetric 
flow rate series and Fig. 16B depicts integrated (or accumulated) injected fuel mass. The first time 
series in both plots relates to simultaneously charging of primary (injector) and secondary coils. 
The second line represents pre-charge scenario. The third curve is the case when charging of SC 
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(AC-wave form of Fig. 1 IB) has been started before the injection (CD-wave form of Fig. 1 IB), 
however, at the moment of 1 .4 ms when the SC-charging was continued the injection has been also 
run. So the overlapping time was 0.6 ms. 

As one can see from instantaneous and integral time series, the fastest opening of the valve 
takes place under shifted charge conditions. The slowest opening is associated with the pre-charge. 
This case also gives lowest level of flow amplitude meaning the lowest speed of the needle at the 
opening moment. A rapid response without any substantial phase delay is associated with the 
simultaneous charge of the SC and the PC. Essentially the same flow amplitude characterizes both 
simultaneous charge and shifted charge. For diesel engines, where the pilot injection and multi-shot 
must be short and produce larger amount of the injected fiiel, shifted charge technique is mostly 
suitable. Simultaneous charge is well applicable to direct injection gasoline engines and also for 
diesel engine at the stage of multi-shot main injection when less stratified fuel spray is desired. 

Some details with respect to each charging scenario at the beginning phases (opening of 
the valve and startup of injection) are shown in Fig. 1 7. There are three plots of instantaneous 
volumetric flow rates along the top row and three plots of integrated (or accumulated) fiiel masses 
along the bottom row. Each of the three correspondents to each of the three different secondary coil 
charging scenarios. The first column reflects data obtained while the SC was simultaneously 
charged with the injector PC (i.e., according to Fig. 1 IB A timing was the same as C timing). The 
second column is related to measurements when the SC was precharged before the injector PC (i.e., 
first was AB of Fig. 1 IB and afterwards started CD, B=C of Fig. 1 IB). The third column shows 
results when the SC charging was shifted with respect to the injector PC operation (i.e., AB and CD 
intervals of Fig. 1 IB were overlapped). 

Under simultaneous charge, the longer the charging time of the SC, the faster opening of the 
valve is observed in instantaneous series as the shift between different series towards the initial 
zero phase. The integrated mass series indicate increased speed of the valve that is seen through the 
slop [g/degree]. The fiiel mean mass rate is characterized by Table 1 below: 



Table 1 - Simultaneous Charge 





M 0.0 ms 


M 1.0 ms 


M 1.5 ms 


M 2.0 ms 


mean mass rate [g/s] 


1.955 


2.07 


2.306 


2.467 


mass per cycle [mg/stroke] 


39.91 


41.4 


46.12 


49.33 



In the case of pre-charge, increasing the charge time results in the same phase of the 
injection startup, but the amplitudes in the instantaneous series and the slops in the integral mass 
series are gradually increasing that says about increased valve speed into the injector. Table 2 
below represents mean mass rates: 
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Table 2 - Pre-Charge 





M 0 ms 


M 1 ms 


M 3 ms 


mean mass rate [g/s] 


0.95 


1.084 


1.122 


mass per cycle [mg/stroke] 


19.01 


21.69 


22.45 



Both effects, the increased amplitude and slopes and more rapid opening resulting in the 
phase shift towards zero phase, which occur under shifted charge technique are shown in the third 
column of Fig. 17. The mean mass rates are in Table 3 below: 



Table 3 - Shifted Charge 





M tau 0 ms shift 0 ms 


M tau 2 ms shift 0.6 ms 


M tau 2 ms shift 0. 1 ms 


mean mass rate [g/s] 


0.443 


0.471 


0.476 


mass per cycle [mg/stroke] 


11.06 


11.77 


11.89 



Application of the SC onto a higher pressure injection system (e.g., over 40 atm of a direct 
injection gasoline system and over 600 atm of a diesel injection system like common rail Bosch) 
results in much more effect on rise time response at the valve opening and fall time response at the 
valve closing. As discussed, for diesel electronically controlled injection system, there may be no 
need to have another SC L2" to rapidly close the valve because fiiel spraying will essentially be cut 
off immediately after first pressure drop. An SC electric circuit consists also of another secondary 
coil L2" shown in Fig. 1 1 A at the position R5. When transistor T2 closes, L2" will produce I- 
Function current in direction opposite to the slowly damping current on the injector primary coil, so 
the resulting magnetic flux will work in parallel with the elastic spring force and results in rapid 
closing of the valve. In another example (which example is intended to be illustrative and not 
restrictive), application of the SC L2" may be important for gasoline and/or direct injection 
gasoline engines where injection pressures are lower than in diesel systems. 

Referring now to the modeling of an electromagnetic actuator according to the present 
invention (e.g., with the second-order non-homogeneous ordinary differential equation (9)), it is 
noted that such electromagnetic actuator ("EMA") may be modeled with an equation different from 
eq.(9): 

x" -haj^x'^ a,,x = a^^ llf\t) - a^, 
by replacing timing components ^'"^ ^ into the right part of the equation to the series of: 



34 



x" + + a^fX = -a^, + 9/ + T + 0, r + .... j ^ 

In this regard, the nature of the added timing derivatives relates to the dynamics of an 
electromagnetic subsystem of a device (or apparatus) to which this particular EMA is applied. The 
coil is ideally represented as an inductor in series with a resistor. In this circuit, the voltage drop 
across the circuit is expressed using the flux linkage dependent on current position of 

the plunger ^ and time phase ' , and coil resistance ^ : 

. d^xj) 
V.^ =r/ + — 

(9.3) 

Circuit current can be expressed as one of the system states by introducing the rate of 
change flux linkage in eq. (93) as: 

dAXxJ) _ dA{x,t) dx ^ dAjxj) di _ ^ 
dt ox dt 01 dt 

The first term ^i^-^'O is determined from magnetic flux ^'"''^^^'^^ 

dx di (94) 

The second term ^2(^>0 is the instantaneous inductance of the coil during transitional 
charge or discharge that can be obtained from dynamic measurements of , ' , , ^/^^ and 
dijdt Because of the parametric nature of such variables, not only the first order of time 
derivatives, but higher orders (second, third, etc.) may be needed to measure and calculate 
regressions to fully construct the right part of eq. (9.1). Of note, from a practical standpoint, 
obtaining an exact analytical solution for the eq. (9.1) may not be possible. However, a numerical 
solution may be found (which implies that on the engineering side it may be essentially impossible 
to have a waveform generator without known input parameters for the electronic circuit). 

Referring now once again to an I-Function, it is noted that such I-Function may take a more 
general form than just the one mode (harmonic) frequency (time) response model of eq. (25): 



J (f \ _^exp(£y22 0 

^fV)"^ (25) 



More particularly, with regard to multiple injection (depicted, for example, in Fig. 19) 
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control over a series of ultra-short injection shots (events) may be utilized for a variety of engine 
operation conditions. Good control of Main! and Main2 may reduce the temperature peaks, and 
hence yield lower amounts of nitric oxides. Pilot shot may yield increased pressure in the engine at 
the end of the compression stroke, thus reducing the start-up time, noise, and smokiness of the 
engine at the warm-up stage as well as increasing the torque at low engine speeds. Pre-M may 
result in reduction of ignition delay that reduces the combustion noise. After-M may provide for 
post oxidizing the exhaust gas and so reduce the amount of particulate matter generated during 
combustion. Post-M is injection of fiiel mainly during the exhaust stroke, thus increasing the 
hydrocarbons HC at the exhaust, and in return, activating and increasing efficiency of the DeNOx 
catalyst. For military vehicles (for example), to increase driving range (fuel efficiency) the first 
three shots, pilot, pre-M and mainl through mainN, may be the most important. The present 
multiple injection driver ("MID") technique may be performed in numerous engineering versions. 
It may be constructed as: (i) a remote electronic driver installed inside a secondary coil; (ii) an 
electronic circuit generating the present I-Function current; and/or (iii) a programmed electric 
current code (e.g., to be incorporated into the main vehicle Electronic Control Unit). 

Accordingly, in connection with a generalized form of the I-Function related to a multi- 
channel MI, each injection shot (event) within an engine cycle may need to be controlled by its own 
channel (e.g., six channels related to the six shot injection sequence of Fig. 19). Each channel may 

have its own time response and phase in order to have flexible control over each 

specific shot (and flexibility in combination of different shots upon the engine run conditions). The 
channels for control of opening and closing the valve may be parallel connections and each channel 
may have a switch controlled by the main Electronic Control Unit that permits a variety of possible 
combinations of the shots. That implies a generalized form for the I-Function as: 



where the primary coil ©21 = 2 ;r K\l Li works in conjunction with a series of secondary coils ©22] = 

2 n R2j/ L2j each of which is switching on and off at its own time phases specified within 
injection cycle. 

Referring now once again to the basic frequency ^1 representing the linear part of the 
complex solution ^ ^» ^ ^ ^2^^' , it is noted that this basic frequency is not, of 

course, solely related to the electrical parameters of the primary coil. Equations (1 1) through (13) 
show what is inside of ^1 , i.e., the normalized parameters in eq. (9) related to the friction, elasticity 
of the spring, gravity and mass associated with all mechanical elements involved into dynamic 
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process (needle, spring, sealing edges, etc.). More specifically, in the eq. (9): 



x'+aj,x' + a„x = a^llf\t)-a^ (9) 



k 



^ - associated mass, 

^lan, . friction coefficient under laminar flow conditions, 

- spring elastic constant, 

- initial elastic force produced by a compressed spring, 
^ - gravity acceleration, 

^« - magnetic field constant, 

- relative permeability, 
^ - number of turns on coil, 

^« - initial spring compression (^*'» ^^), 

- current amplitude, 

^fr ' ^e/ 9 ^mag ' ^sys . transfomi coefficients. 

So, , related to the ^^^^^ " ^ »^^' solution of the first, linear part, represents all 

mechanic, hydraulic and elastic elements of the system while ^J^'^^^'^'^ , related to 
the ^2^^^ = ' solution of the second, non-linear part, represents the parameters of the system 
under impact of magnetic flux. 

Referring now to time-dependent action (e.g., movement of various physical elements) and 
frequency-dependent action (e.g., movement of various physical elements) of the electromagnetic 
actuator (e.g., dependent upon the resistance R2 and the inductance L2), it is noted that a 
generalized impulse balance has been identified in eq. (1) as: 

Now, consider the moment at which magnetic force becomes over all others involved in the 
process. From this moment the equation can be simplified to: 

dU „ 

dt (l.l) 

To derive a relationship between the velocity ^ of lifting armature (or valve, or needle, or 
associated mass in general) and transient current (I-Function), one needs to make energy balance on 
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electromagnetic and electric parts . That can be performed in terms of power release: 

Wmech=Kn, (1.2) 

Mechanical power is work ^ over time , so using impulse, it can be expressed as: 

dA Fdx dU d(mU') 

dt dt dt 2dt (J 3^ 

The voltage over the coil is dependent on current derivative: 

dt (1.4) 
Electromagnetic power is related to instantaneous voltage and current: 

dt 2 dt 

In the case of balanced energy transfer, the relationship between lifting (pulling in/out, 
pushing in/out the armature) velocity and current time series becomes linear: 



^ (1.6) 



This equation implies that in order to get control on rapidness of the primary solenoid 
(injector solenoid) with known inductance ^> and associated mass ^ , the speed and transient shape 
of the lifting is directly related to the current time series. The acceleration ^ (or force is 
proportional to the first order current derivation: 



dt\m 



(1.7) 



Eqs. (1.6) and(1.7) are very important for both injectors and electromagnetic air valvetrains 
to control speed-acceleration control during opening and closing the valve. In the case of fiiel 
injectors both the opening and the closing events must be rapid in order to make stability (e.g., 
gasoline injectors) and/or multiple injection (e.g.,diesel injectors) possible. In the case of air intake 
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valve, the rapidness (maximum speed and acceleration) are important at the opening of the valve, 
however, by closing the valve at the end of armature movement, the speed and acceleration must be 
close to zero (problem of durability). 

In this regard, the diagrams in Fig. 10 represent lifting velocity (top diagram) and 
acceleration/deceleration (bottom diagram) for three different ratios between primary and 
secondary coils in arbitrary units. 

For the primary coil, the angular frequency ^ ^1 is represented as series of 40, 1 5 
and 5 units. For the secondary coil, its frequency ^22 = 2;r I \^ represented as series of 20, 10 
and 5 units (always slower). The higher the ratio ^^22 , the higher the rapidness in both terms of 
velocity and acceleration. 

The time phase where (*^'^^^)22of the secondary coil becomes the minimum is a time phase 

when the transfer of energy from secondary solenoid to primary solenoid should be ended. This 

time ^22 has to be equal or proportional to the time response of the whole dynamic system ^^^""^ as 

it sketched in Fig. 14, which is determined by injection combustion conditions. For example (which 

example is intended to be illustrative and not restrictive), to make possible diesel multiple injection, 

the dynamic rise/fall time should be not longer than about 200 u s. To achieve that, under this 

example, the electromagnetic actuator (primary coil) must react in about 100 w s. The factor of 
tit < 1 

22 dymm verified experimentally (e.g., using the instantaneous fuel flow rate technique 

discussed herein and/or high speed visualization of the fuel spray). So, the final setup of ^22 is an 
iterative process starting from a lower ratio of ^21 ^^22 and incrementing it until the value of 
will be within a given range. 

Referring now to how the time-dependent action and/or frequency dependent action of the 
electromagnetic actuator may be determined (e.g., calculated, measured), it is noted that one 
example algorithm (which example is intended to be illustrative and not restrictive) is described 
below. More particularly, this example algorithm of the determination of time response (J^^^Jn ), 
frequency (^22)^ and coil (^2»^2) is as follows: 



• Cycle # 1 - Construction of the secondary coil driver ("SCD"). 

1. Upon engine model, injection system model, fuel load map at different engine runs 
(speed versus torque-horsepower) timing strategy, exhaust emission requirements and 
electric configuration (ECU injection timing control, RL characteristics of the injector 
solenoid, voltage applied, etc.), the first injection pattern is designed such as shown in 
Fig. 19, specifically: 
o Number of shots, 
o Duration of shots 
o Rise/fall times. 

o Dwell intervals between the shots 

o Fuel amount per shot (amplitude profiles). 
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o Tolerance range for time phases and amplitudes (fuel amounts). 
One would understand that Fig. 19 could hypothetically form the basis of a 
corresponding curve having time on the x-axis (in arbitrary units) and current on the y- 
axis (in arbitrary units). 

2. Determination of the using instantaneous fiiel flow-metering technique, 

3. Limitation of ^ 

4. Determination of ^22 by doing numerous iterations to get curves of the I-Function at 
^22 within given tolerances (time and amplitude). Of note, the iterations generate curves 
which can be compared to values of Fig, 19; the curve closest to the one capable of 
producing Fig. 19 indicates the value of ^22 

5. Knowing lifting velocity ^ " ^'^^ ^^^""'"and ^"^^p^^ calculate ^2 using eq. (1 .6). 

6. Calculation of ^2 = ^22 4 /(2^) . 

7. Construction of secondary coil driver (as physical unit or electric circuit or programmed 
I-Function code) with variable A . 

• Cycle # 2 - Testing of Multiple Injection with Applied SCD. 

1 . Test of injection pattern under various injection cycles (frequency, number of shots, shot 
duration, dwell intervals) to see output dynamic characteristics using instantaneous flow 
metering technique. 

2. Repeat Cycle # 1 to achieve required rapidness and stability. 

3. Test injection system in long run (-100,000 cycles) to validate durability. 

• Cycle # 3 - Engine Test. 

1 . Install injectors in engine equipped with SCD between injection timing driver and 
injectors. 

2. Test engine performance (power and torque release) to achieve max fuel efficiency at 
the required torque output using a dynamometer test cell. 

3. Test engine exhaust emissions, 

4. If necessary, repeat Cycle # 2 to change injection pattern as required. 

5. Test engine in long transient and steady state run. 

• Cycle # 4 - Road (Drivability) Test (Extended). 

1 . Install injectors on a vehicle with the same injection system, which has been tested 
during Cycle # 3. 

2. Measurements of the fuel consumption (continuously) and exhaust emission (selected 
test) at different driving and climatic conditions. 

3. If necessary, repeat Cycle # 2 to change injection timing/phasing strategy to minimize 
fuel consumption and exhaust emission. 



40 



Regarding Cycle #1 above, it is noted that in this example the phasing of the I-Function 
itself and its peaks are related to Fig. 19 in that Fig. 19 represents the injection mapping target upon 
certain engine demand (i.e., regime). 

Further regarding Cycle #1 above, it is noted that in this example ^'^^'"is determined on the 
basis of measured time series of instantaneous flow rate along with velocity, pressure gradient and 
integrated mass series. To determine this time factor one can use either flow rate or pressure 
gradient time series. In the first there is a dynamic rise sharp slope which is ended by a zigzag-type 
peak. This peak says that the valve is opened, the injection has actually occurred and the break-up 
point (transfer of the liquid jet into droplets) has taken place. The angle of this slope represents the 
speed of this dynamic process, i.e., how fast the whole system (mechanics, hydraulics and inertia of 
all associated masses) has reacted after a given electric wave form onto the primary coil (injector). 
In the series of pressure gradient this factor is determined by a rapid spike-like change of pressure 
gradient from negative (acceleration of the flow) to a positive derivative. 

Further regarding Cycle #1 above, it is noted that in this example lift of the injector valve is 
a design property which is essentially a fixed parameter. For instance, in direct injection gasoline 
engines it is typically about 50 to 90 micrometer, in normal gasoline injectors it is typically up to 
about 300 micrometer, and in diesel injectors it is typically between 100 up to 500 micrometer. In 
other words, //// is a given parameter which represents a gap between a sealing position and a 
pushing upward/downward stop position. 

Referring now to another embodiment of the present invention regarding an application 
related to controllable high-pressure fuel injection in diesel and direct injection gasoline engines by 
means of stable multiple ultra-short injection events using a secondary coil driver (SCD), attention 
is directed to Fig, 18 (such multiple injection under stable timing and amount controlled by SCD 
provides a cascade-like fuel spray and flame structures with a more widely spread surface for the 
compressed air, as depicted in Fig. 13, for example). Of note, an important element in such an 
injection technique is the events' (shots') timing that may need to maintain a core flame to prevent 
a quenching effect. Thus, the final spray structure may have the appearance of a flipped-cascade 
Christmas Tree in which only the jet and premixed zones are fully developed without the 
appearance of the rich zone. 

In this regard, the combustion process in reciprocating internal combustion engines is a 
complex dynamic phenomenon including ftiel injection, air intake, air-fuel mixing flow, chemical 
and thermodynamic kinetics, mixture burning, and exhaust of combusted gas with pollutants. This 
dynamic process has different time scales in terms of the engine in-cylinder kit reciprocation, fliel 
injection, chemically inter-reacting species kinetics, ftiel spray and flame formations. All these 
timing scales become extremely important in high-pressure injection engines such as diesel and 
direct injection gasoline engines. 

More particularly, the reciprocating cycle fits an order of a few tens of millisecond (--lO'^ 
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sec). Injection lag is about a few hundreds of microseconds (-10"^ sec), and injection duration has a 
few milliseconds (-10'"^ sec) in gasoline engines. In diesel engines injection lag and injection 
duration are shorter, --10'^ sec and -10"^ sec, respectively. In local flame domains, the ignition lag 
and premixed flame and rapid oxidation (combustion) in diesel engines have an order of magnitude 
of a few tens of microseconds (---10"^ sec). In gasoline engines these factors become a few hundreds 
of microseconds (-lO""* sec). Typically, in diesel engines all processes are more rapid having one or 
two orders shorter duration. 

An important conclusion is that injection shot Atsh and dwell duration Atdw may have to be 
directly related to the early stages of diesel combustion, i.e., in the manner of timing of injection 
dynamics and chemical kinetics (in the case of single shot per cycle, the sequence may begin 
shortly after the start of fiiel injection and may continue through the premixed bum and into the 
start of quasi-steady combustion). 

The time between the start of injection and the premixed bum may be about a few hundred 
microsecond (-10"^ sec). If, at that moment injection stops, the premixed zone may start to be 
developed in that space and completely burned as a regular premixed reacting substance. This 
factor may determine dwell interval to be close to •-'100 u sec in order to exclude in the combustion 
process a further development of a fuel-rich zone. 

The injection ultra-short shot duration may be determined by the time limit needed to get 
the injection of about -1 u sec started, i.e., by injection lag. Depending on the fuel amount demand, 
the production factor may be varied, for example, from about 10 to 30, meaning that shot duration 
in this example may be about --10 to 30 u sec. 

In another example (which example is intended to be illustrative and not restrictive), the 
exact set up Atsh and Atdw for a particular type of engine and injection system may be dependent 
on: 

1 . Fuel properties such as density, kinematic viscosity, surface tension, boiling temperature, 
specific heat and/or compressibility factor. 

2. Injection pressure fluctuations. 

3. Nozzle geometry. 

4. Compression ratio. 

5. Partial fuel load per cycle. 

Thus, the need may arise to test a fuel injection system and engine at different loads and 
speeds to tune the SCD for the final setup of Atsh and Atdw at different mapping conditions. To 
make the SCD work in conjunction with a certain type of engine and injection configuration, it may 
be necessary to proceed with the following example subsequences (which example is intended to be 
illustrative and not restrictive); 



1 . Analysis of high-pressure injection dynamics (an OEM's original injection system) by means 
of the instantaneous flow rate measurements indicating exact positioning of breakup peaks and 
ICCD (Intensified Charge Coupled Device) high-speed fiiel spray visualization in order to 
prove spray structure both in terms of liquid (fuel jets and droplets) and gaseous (evaporated 
fuel) phases. 

2. Design, simulation and construction of a secondary coil driver (SCD) applicable to a 
production injection system. 

3. Experimental verification of rapidly controlled multiple injection by means of flow rate and 
fuel spray dynamic measurements as in Step # 1. 

4. Experimental verification of in-cylinder diesel fuel mixture without and with applied SCD. 

5. Tuning of engine performance and emissions in one-cylinder engine model without and with \ 
applied SCD. 

6. Tuning of the OEM's engine performance and emissions in a production model without and 
with applied SCD under the tuned discharge method. All engine torque-speed operational 
diagrams may need to be mapped. 

7. Design, construction and test of industrial SCD on-board prototype either in the form of SCD 
or electric circuit or encoded I-Function current electronics. 



Referring now to Fig. 19, certain injection events associated with an example of the present 
invention (wherein the injection events are identified with reference to certain combustion effects 
and engine run/injection strategies) are depicted. More particularly: 

• With reference to certain combustion effects 

o M 1 M2 . . . -reduces T-peaks (NOx), fuel consumption 

o After-M - provides post oxidizing exhaust gas (PM) 

o Post-M - increases HC in exhaust (DeNOx catalysts) 

o Pre-M - reduces ignition delay (noise) 

o Pilot - increase P in cylinder (start-up, noise/smokiness at warm-up, torque at low 
speed) 

• With reference to engine run/injection strategies 

o Engine start-up/warm-up: Pilot - Pre - Mainl 

o Texhaust < Tcataiyst : Pre - Mainl - After 

o DeNOx TEC: Pre - Main 1 - Main2 - After - Post 

o High TEC: Pre - Mainl - Main2 - After 

o High Torque, low speed: Pilot - Pre - Mainl 

o Medium/high speed & load: Pre - Mainl - Main2 

o Maximum-power conditions: Pilot - Main 1 



43 



Reference will now be made to an example (which example is intended to be illustrative 
and not restrictive) of certain engineering calculations to design a secondary coil and 
coding electric current to be applied to an injector (e.g., a Bosch common rail injector). Of note, 
this example is aimed at a simple demonstration of what needs to be known, calculated, coded, and 
transferred to a primary solenoid actuation variety of device. This particular example is directly 
associated with a production Bosch common rail injection system (CRIS). A commercially 
available L/C Meter IIB in the u H range has been used to measure inductance of each of four 
injectors installed on the CRIS. An HP/Agilent 33 120A 15 MHz Function/ Arbitrary Wave 
Generator along with HP3481 lA BenchLink Software are applied for output signal coding of the 
voltage/current time series. And HP Infinium 500 MHz IGsa/s Oscilloscope has carried out 
verification of quality and time phases of the output control signal fed to the CRIS injectors. 

In summary, the algorithm steps described below can be divided into three basic stages: 

1. One needs to measure electric properties of the injector, such as inductance L and 

resistance R , to evaluate time/frequency response. That permits a calculation on energy 

transfer per each peak, spike or other fraction of the injection timely controlled electric 

current/voltage. Now, at a given factor of energy transform, it becomes possible to 
R f 

calculate ' - parameters of the secondary coil (SC) which must generate a transient 
current to make rapid opening and closing of the valve. 

2. Now, one needs to proceed on I-Function current as a time series and determine what 
time phase (charging time) is most applicable for rapid, stable control over actuator. For 
instance, with regard to gasoline injectors or diesel injectors with electronically 
controlled hydraulic valve, at the valve opening stage the part of time series may range 
from the beginning until a phase where the I-Function current has the maximum because 

instantaneous velocity of armature is proportional to instantaneous current " ~ ' VV^ ^ 
In the case of an air intake valve it may be necessary to have the time series until the 
first current derivative becomes almost zeroed. This is due to proportionality between 
instantaneous acceleration (force) and current derivative ^ ~ (di^dt)"^ yfZ/m ^j^^ 
will be incorporated into an injection system as firmware, at this stage the algorithm 
may switch to an electric fabrication of the SC driver and tuning it in terms of discharge 
mode (described above). If the SC is to be implemented as a code, the process continues 
to the third stage (below) 

3. The obtained I-Function current time series may be fitted to a standard library function 

available in software to an arbitrary (ARB) wave generator. Now, after matching the 

R f 

derived I-Function upon ' -characteristics of primary and secondary coils and the 
librarian one, the setup of mathematical parameters becomes available to construct 
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different transient phases of the injection cycle including individual injection shots and 
their u s time fractions. Finally, the constructed current code may be transferred into the 
given ARB-generator that controls the injection profile. This procedure may need to be 
repeated a number of times to cover an OEM's injection map. Afterwards, it is possible 
to transfer the entire SC driven injection map into a processor that is incorporated in the 
vehicle ECU. Depending on driving and engine run condition, the ECU may call either 
the OEM's or the ARB injection control current codes related to a particular injection 
event in each injector. 

Referring now to aspects of the detailed algorithm outlined in the 3 stages above: 

1 . OEM's Injection Map It may be critical to know the exact technical data regarding the 
OEM injection system, injector operation, and current/voltage trace applied onto the 
actuator may be required. The solenoid valve (triggering element) may control a valve 
ball and at the stage of its pulling in (energized solenoid) the bleed orifice may be 
opened (and a pressure difference between the feed passage to the nozzle and valve 
control chamber causes upward lift of the nozzle needle ~ resulting in injection event). 
The energizing time of this solenoid varies (e.g., from 1 to 2 ms) with a peak puiling-in 
current of, for example, 18A and holding current of, for example, 12 A. The rise time 
and fall time are varied (e.g., from 80 to 100 w s). During the holding stage the current 
oscillates (e.g., with amplitude of 0.57A and periodicity of 0.1-0.2 ms). A typical 
current trace applied to the Bosch CRIS injector is illustrated in the left plot of Fig. 25. 

2. Actual Injector Solenoid RL Data Resistance R was measured using a multimeter. 
Inductance^ has been obtained using L/C Metter IIB that has a wide range of 

^ sensitivity from nH, u H, mH to H. Zero mode has been constantly applied to subtract 
the stray inductance which was initially about 1.8-2.2 u H due to measurement wiring 
and after Zero mode was oscillating at 0.007 u H due to the wire loop configuration and 
temperature dependence of resistance during measurements. Referring now to Fig. 20, 
the RL data are shown along with time and frequency response characteristics of the 
injector (primary) coils. According to both measurements, and the left plot of Fig. 25, 
the rapidness of different solenoids (rise-fall time) at the opening and closing the valve 
is varied from 146 to 212 w H (resulting in a frequency response of 4.72 to 6.85 kHz, 
respectively). In two columns of Fig. 20 the power ^ )M ^ fluxed into the 

primary solenoid during energized state is calculated using measured inductance ^ , 
pulling-in peak ^^'"^ - 18^ holding "^^"^ current, time response and holding 

duration respectively to peak and holding stages. As indicated, p^"'' varies from 64.8 
p 

to 72.9 W and '^^'^ =4.7-6.1 W for various injectors. These power (energetic) values may 
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be limited by construction of coil, i.e., its inductance ^ and currents ^^''* , ^'"''^ upon 
dynamic time response. 

3. Targeted Power and Time Response Conversion Ratios. SC ^ " ^^^^ To make the 
solenoid more rapid that results in stable ultra-shot injection needed for controllable 
multiple injection it may be necessary to have additional energy that will be released 
very rapidly may be required. In the Bosch CRIS the electromagnetic actuator 
(solenoid) controls opening and closing the valve. The distance between the high 
pressure inlet into the injector from CRIS to nozzle needle chamber is 0.11 m, the sound 
speed under 1600 bar is -1700 m/s, so the time of pressure propagation is about 65 u s. 
That implies a magnitude of time fraction that must be comparable with minimal 
rise/fall time of the actuator and very stable (repeatable) from cycle-to-cycle. The 
secondary coil does produce a quick additional power to speed the rise/fall phases. In 

the right part of Fig. 20 the calculation of ^ -parameters are reflected. The first input is 

E E 
the power ratio between p^""^^ of the primary coil and ^''"^^ of the secondary coil 

E = FE 

peak! peak! ^ ^hg^c factor F is varicd between 1.5 to 4.0 depending on the actuator 
type and its application. In this particular example, its maximized because for multiple 
injection in diesel injection with "light" inductance (high response time) the effect of 
rapidness is associated with high power ratio input ^ = 4.0 , That permits calculation of 

inductance of the secondary coil as p^"^^ ''^^^^ '^'^ ^ . Reversibly, the secondary 

coil has slower time response z'^^^^- peaki ^ ^^ere 2.0 < ^ < 5.0 Once again, because 
multiple injection requires a quick control over both injection shot and dwell interval 
between injection shots, factor ^ - 2.0 is minimized. That results in resistance value 
R = L IT 

2 2 peak! ^j^c SC drivcr is desired to be made as a physical electronic circuit, 
R^L^ - data enough to design and construct as described above. If the I-Function 
current must be driven as wave-shaped code, it may be necessary to proceed on to the 
next four steps. 

4. Construction of I-Function Having frequency responses of both primary and secondary 
coils, one can construct an I-Function current timely trace in normalized to unit form as; 

\2TcRylL^]t 43.0? 

i{t) = /exp^^P^f^^^^^^ = lexp'^^P^^ 



Such an I-Function current trace and its first derivative are shown in Fig. 21. Because 
^1 L data are in kHz, the time scale is in ms. The maximum current peak corresponds to 
0.047 ms which is related to the maximum velocity of the primary solenoid armature. 
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That time duration is a time '^^'^^^ that should be given to the secondary coil to be 
charged before transferring the energy to the primary coil. 

5. Fitting to Library Standard Waveform Waveform generator hardware can reproduce a 
variety of current traces calling so-called standard waveshapes and their combinations. 
That moves the algorithm to the next step, which is translation of an I-Function current 
into available library functions and time into number of points within the cycle. For 
example, in HP 33120A software one cycle is equal to 16000 points (pts). For the rise 
and fall I-Function current most fitting shapes are rise ^(^ " ^ ) and fall 

exponential functions. In normalized to unit form, the amplitude ^ is equal 1 . So, 
the damping fact ^ maybe be derived from comparison with I-flinctions at rise and fall 
fractions: 

1-exp''*-" =l-exp"*" 
exp'"*"" =exp'''" 

where K, Q, and n are determined during the fitting process (the resuh of which is 
shown in Fig. 22). 

In this example, one has the following equations: 

• 0.175m. 2.25 
0.047m5 



02132!£iii = 9.60 
0.047mj 

6. Targeted Multiple Injection Map and Time Scaling Fig. 23 indicates translation of 
camshaft angular positions of various phases during an injection cycle. In this example, 
the engine speed is 400 RPM for four stroke cycle (f=33.33 Hz). The main injection is 
set at 180° (top death center TDC). Before TDC at -20° starts pilot injection. Both shots 
have duration of 600 u s. Dwell interval is 1275 u s. All phases are calculated in degree, 
u s and pts. 
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7. Construction of Special Waveform Each phase may be coded. Fig. 24 illustrates the two 
shot injection per cycle calculated in previous step 6. As shown, each shot is divided 
into 5 phases and translated into absolute and arbitrary coordinates of time and 
voltage/current amplitude. The resulting output signal is shown in the right plot of Fig. 
25. 

In another embodiment of the present invention Angular frequency ^21 = 2;r7?, /Z, ^Yad/s]; 
Frequency -^21 ^ A [Hz]; Time response (rise) ^21 = A /^i [s or ms or u s]; Angular frequency 
0)22 = 2^/?2 [rad/s]; Frequency -^^^ [Hz]; and Time response (rise) ^22 =^2/^2 [s or ms 

or ms]. 

In another embodiment the present invention provides for application of I-Function ultra- 
short transient magnetic flux cutting transient inertia in wave form diagrams of solenoid-valve 
needle stroke (or more generally, coil-plunger stroke) that results in rapid dynamic of 
force-stroke response (solenoid performance). 

In another embodiment the present invention provides for theoretical solution(s), actuation 
technique(s), engineering realization(s) and/or experimental method(s) related to rapidly operated 
injection. 

In another embodiment the present invention provides an exact analytical generalized 
solution to a second-order non-homogeneous ordinary differential equation describing complex 
dynamics in a primary solenoid including magnetic flux, elastic force, gravity and friction. Of note, 
this solution indicates that spectrum characteristics (frequency and/or time response) are fully 
dependent on time-dependent transient current applied at the opening and closing of the injector or 
any other like actuator. This current can be generated from an outside source (outer from 
primary solenoid). 

In another embodiment the present invention provides an "I-Function" which satisfies a 
frequency and/or time response relationship between a remote secondary coil and a primary coil in 
terms of resistance to inductance ratios. Of note, the strongly exponential I-Function has unique 
features that help determine main criteria to construct secondary coil and/or a current electric 
circuit to the drive primary solenoid in an injector or an actuator. 

In another embodiment the present invention provides inductive pre-and post- secondary 
inductor circuits for a fuel injection system or any other like actuator in order to control both 
rising and falling time response at the opening and closing of injector valve (or in more general 
application, the plunger opening and closing dynamics related to an electromagnetic actuator). In 
one example (which example is intended to be illustrative and not restrictive), this circuit may be 
flexibly constructed for wide application range by changing nominal characteristics of different 
circuit components with respect to a particular application case on the basis of primary solenoid 
characteristics and/or time response limits needed for injector or actuator rapid operation in a real 
environment. 
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In another embodiment the present invention provides at least two different secondary coil- 
charging techniques (referred to in the present application as simultaneous charge and pre-charge). 
Of note, these different charging scenarios indicate that transient I-Function current can be shaped 
in different ways in order to manage its aptitude-time-spike wave forms for different actuators. In 
another embodiment the shifted charge technique, which is combination of the first two scenarios, 
is also realized. 

In another embodiment the present invention provides instantaneous fuel flow rate 
measurements applied to indicate that the remote secondary coil technique not only generates rapid 
electric I-Function current, but also results in rapid transient dynamics in the instantaneous flow. 
Such instantaneous fuel flow rate measurements support certain theoretical and engineering 
conclusions discussed above. 

In another embodiment the present invention provides that the I-Function may be generated 
from the secondary coil driver without physical usage of the coil. That is, the I-Function relates to 
a current to be applied onto a primary solenoid in an actuator. In another embodiment an I-Function 
current generator may be utilized knowing basic parameters of primary solenoid. Such a current 
generator (or driver) may produce current to be applied in the form of a time-series coded 
waveform (e.g., from a resistor to which time-dependent voltage is applied). 

In another embodiment the present invention provides that the I-Function may be directly 
coded (e.g., as a binary code into a chip installed into an Electronic Control Unit of a vehicle). 

In another embodiment the present invention provides that the I-function may be coded as 
software. In another example (which example is intended to be illustrative and not restrictive), such 
software may be transmitted (e.g., through the Internet) to a solenoid to operate a remote actuator 
within given time limits of its opening and closing stages. 

In another embodiment the present invention provides that the I-Function control technique 
may permit improvement in time response characteristics of existing devices in industries where 
timing is important for the whole dynamic process. In one example (which example is intended to 
be illustrative and not restrictive), application may be to diesel engines (to permit control of 
multi-shot injection as a series of ultra-short pilot injection and multi-shot injections within main 
injection as well as to control dwell interval between injection shots in order to get complete 
combustion and ultimately decrease fuel consumption and emission of particulate matter and 
nitrogen oxides (i.e., high injection repetition rate controller)). 

In another embodiment the present invention provides for increasing vehicle fuel efficiency 
(e.g., diesel fuel efficiency) and/or driving range of vehicles equipped with either common rail or 
unit injector or unit pump or distribution injection pump systems. 

In another embodiment the present invention provides for a multiple injection driver (MID) 
to implement controllable and timely repeatable multiple injection. 

In another embodiment the present invention provides for a controllable injection phase 
shift (e.g., advanced and/or retarded), in order to get efficient and complete combustion and 
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heat/pressure release. 

In another embodiment the present invention provides for the utilization of existing serial 
electromagnetic actuators mostly constructed by using a single coil assembly. Analysis and 
realization of their rapid switch on/off operation essentially without transient delays are carried out 
with reference to Figs. 6A-6D and 7A-7D, for example. More particularly, one or more of the 
following may be utilized: 

• Analysis of the transient mechanic and electromagnetic dynamics which typically occur 
during an electromechanical actuator operation (with focus on the start/end transitions). 
This part considers general theoretical analysis through representation of an exponential 
type time dependent solution obtained under the gravity, magnetic, elastic and friction 
forces applied onto the injection valve. 

• Introduction of an I-Function, which is generated by a remote secondary coil in the from of 
a rapid transient induction current to be applied onto primary solenoid. 

• Engineer an electric circuit to realize SC-technique with respect to internal combustion fuel 
rapidly operated injectors. 

• Realize a program which computes charging (energizing) time of the SC under defined 
properties of the PC. 

• Experimental verification, including electrical measurements and measurements of 
instantaneous fiiel flow rates, indicating simultaneously the complex dynamics of 
electromagnetic, hydraulic, mechanical and frictional factors contributed in final time 
response of the injector. 



EXAMPLES OF PERFORMANCE TESTS AND QUANTIFICATION ACCORDING TO 
EMBODIMENTS OF THE PRESENT INVENTION 

I. Performance Evaluation of a Multi-Burst Rapidly Operating Secondary Actuator applied 
to Diesel Injection System 

INTRODUCTION 
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The following now refers to a performance evaluation of a multi-burst rapidly operating 
secondary actuator according to an embodiment of the present invention as applied to a diesel 
injection system. This embodiment of the ROSA is aimed at further improvement of diesel fuel 
efficiency and exhaust emissions. In this regard, the inventor has conducted tests of ROSA 
aimed at providing controllable and repeatable multiple injection events, particularly in common 
rail injection systems ("CRIS"). Currently, fuel system suppliers are typically resorting to 
piezoelectric switches and other costly electric and electronic control units to provide the multi- 
firing effect in CRIS. ROSA generates a special current, which is applied onto the primary 
solenoid of the injector to control its transient fast response. An injection test cell has been 
constructed for this performance evaluation. Two test setups were available for both diesel spray 
visualization and instantaneous fiiel flow rate measurements. Up to six shots per cycle were 
implemented under injection pressures from 1200 to 1800 bar. The injection repetition rate was 
equal to a four-stroke engine speed of 1200-3600 rpm. A high-speed digital camera was used to 
have accurate quantitative data regarding diesel spray rapid dynamics. An argon laser 
illuminated the spray field. Processed data were obtained for liquid spray tip velocity, injection 
shots duration, and their delay with regard to electric signal setup. The stability of phasing lies 
within 50 |as. The shortest injection shot duration is 74 (as, maximum variability of short duration 
is 50 \xs. An advantage of ROSA is very stable phasing, dwelling and duration of multiple 
injection shots proved from cycle-to-cycle analysis. The ROSA technique also has a number of 
other unique applications including Electronic Unit Injector (EUI) and Hydraulic Electronic Unit 
Injector (HEUI) and variable air intake valve actuators. 

Recently, it was shown that multiple injection technique, applied to different diesel injection 
systems, has tremendous practical potentials to improve diesel combustion and aftertreatment 
processes in variety of engine performance characteristics, including fuel consumption, 
emissions of soot/NOx and noise. There are numerous strategies in the split of single main 
injection into a series of sequential events, namely called Pilot, Pre-Main, Main-1 and Main-2, 
After-Main and Post injection event or shot. They can be summarized as is illustrated in Fig. 28 
for six shot injection with arbitrary referenced cam phases within the injection cycle. For 
instance, good control of the Main injection(s) reduces the temperature peaks, and hence yields 
lower amounts of NOx. The Pilot shot yields increased pressure in the engine during the 
compression stroke, thus reducing the start-up time, noise, and smoke level of the engine at the 
warm-up stage, as well as increasing the torque at low engine speeds. The Pre-Main injection 
event results in a reduction of ignition delay thereby reducing combustion noise. The After-Main 
shot provides oxidization of the exhaust gas, which reduces the amount of particulate matter 
generated during combustion. The Post injection occurs during the exhaust stroke, thus 
increasing the hydrocarbons HC at the exhaust, which increases the efficiency of the DeNOx 
catalyst. Most of the multiple injection studies are directly related to the CRIS type injection 
systems. Fewer studies are highlighted with regard to EUI and HEUI, mostly applied for heavy- 
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duty diesel engines. 

To make multiple injection systems widely practical in automotive industries, it is necessary 
to provide very stable timing associated with four factors. The first is phasing of injection shots, 
the start of injection events. The second is injection duration of each event. The third is dwell 
interval between shots, especially related to Pre-Main, Main-1 and Main-2. And the fourth is 
delay factor dealing with the time needed for pressure propagation along the high-pressure pass 
from a pressure accumulation or generation source to an injector control valve as well as for 
pressure recovery. All these timing factors become very critical in the following cases: (i) 
increased number of shots, e.g., up to six; (ii) shorten dwells, e.g., down to 200 ^s; (iii) enlarged 
dynamic (max/min) range of injection fuel flow rates for different shots, e.g., --100 mg per Main 
and ---0.1 mg per Pre-Main; (iv) uncontrolled fiiel pressure oscillatory frequency (-10-100 Hz) 
that can be in resonance with some multiple injection harmonics. These harmonics are widely 
varied from a few Hz to a few kHz. 

As can be seen from various engineering conceptual designs of injectors and injection 
systems applied for multi-firing, there are one or two valves that control fiiel pressure 
distribution between control and accumulation volumes associated with spill and needle valves, 
respectively. In older injector generations such as 1^^ generation of CRIS an electromagnetic 
actuator controls a spill valve, which is hydraulically connected to a high-pressure line fed 
directly to the common rail (almost constant high-pressure source). While triggering the injector 
spill valve by energizing a solenoid type actuator, the pressure in control volume drops down 
below the pressure in accumulation volume. When the pressure difference applied on the sealing 
area of the injector needle overcomes the needle spring force, the injection starts. So, the 
actuation of injection in such solenoid type electronically controlled diesel injectors is a one- 
stage process. In some systems, where the piezoelectric actuator or second actuator (for instance, 
two-actuators EUI) hydraulically coupled to the needle valve in relatively closer position to the 
needle spring, the timing control on fiiel pressure propagation to the accumulation volume can 
be flexible split into two stages. 

At the first stage, the spill valve controls pressurization of the entire high-pressure gallery of 
injector by a common rail in CRIS or a pumping plunger in EUI or HEUI. Then, at the second 
stage, the needle valve controls the injection process itself Practical implementation of new 
muhiple injection techniques is quite costly and cannot be applied to the series of existing 
electronically controlled diesel injectors. 

Only a few studies related to the timing stability of multiple injection are currently 
available. For example, cycle-to-cycle variability in injection characteristics was observed and 
explained by cyclic pressure deviation up to 22% in the common rail. Different timing strategies 
for the split of main injection into Pilot, Main and After with shifted phase and duration are 
studied, but only constant delay of the actual injection relative to the electric trigger signal of 
about 300 ^is is outlined as a factor of stability. There is also a little data related to well- 
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quantified amounts of fiiel injected per each shot. Regarding a production multiple injection 
system, up to a 5-shot system with 400 dwell between Pre-Main and Main events and 
minimal injection fiiel amount of 1 mmVshot with controllable variability of 0.5 mm^ was 
mentioned in 2003. 

The present inventor has developed a novel technique for a variety of applications related to 
the rapid acceleration and deceleration of a plunger into an armature, where the high timing 
stability is crucial for a specific process. With regard to automotive applications, primarily 
applied on any electronically controlled fuel injectors and variable air intake valves, this 
technique is based on a rapidly operating electromagnetic secondary actuator (ROSA) triggering 
the pressure control valve solenoid installed onto/into the injector. Physically, ROSA generates a 
specially shaped current called I-Function current, which is transferred onto the primary solenoid 
of the injector. This current controls the rise and fall transient response of the primary solenoid 
that results in controllable rapid and stable opening and closing of the injector valve. 

The ROSA technique can be performed in numerous engineering versions including (i) a 
remote secondary coil (for medium- and heavy- load solenoids of injectors and air intake 
variable valves for diesel engines), (ii) an electronic circuit (for lower load devices such as 
gasoline injectors), and (iii) a coded current profile incorporated into vehicle ECUs/EDUs. In 
this particular project, an in-coded version of ROSA was constructed and applied to a first 
generation Bosch type CRIS designed only for single shot injection with min/max energizing 
duration of 1-2 ms respectively. The main objective of this study was a quantitative validation 
of ROSA multiple injection control by means of a high-speed visualization of the diesel spray. 
In this case, the operation of the entire injection system results in a spray dynamics out the 
injector as shown in Fig. 29, Accurate temporal and spatial recording of the spray sequences 
provides detailed information about fast transitions occurring during high-pressure injection. The 
temporal resolution must be close to a few tens of microseconds to observe primary break-up 
transition, jet tip supersonic velocity and all injection timing characteristics needed for the 
required validation. 

Details of the performance evaluation are as described below: 

ROSA-CRIS EXPERIMENTAL SETUP 

General Configuration 

Initially, the utilized CRIS was not equipped with a production electronic control unit 
(ECU). A Kistler 4067A2000 piezoresistive high pressure sensor along with a 461 SAO amplifier 
measured pressure in the common rail, which was without a pressure limit switch to control the 
CRIS spill valve solenoid. Fig. 30 illustrates the technical stages that were carried out in order to 
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construct an integrated test cell. Four subsystems, i.e., (i) a high-pressure (HP) hydraulics unit, 
(ii) a ROSA based electronic injection driving unit (EDU), (iii) a volt-to-amp converter, and (iv) 
a high-speed visualization channel have been constructed and incorporated into the test cell. The 
interconnections between all subsystems are shown in Fig. 31 along with the specifications of 
equipment used. The system allows very flexible and fully controllable setups of input and 
output data using two PCs. 

High Pressure Hydraulics 

The HP hydraulics unit is composed of a 40-liter fuel tank, a low-pressure pump with a fuel 
filter, a high pressure 5 Dm- filter, an electric motor which motorizes a high pressure pump 
connected directly to the CRIS. An additional electric controller was used on the motor to have a 
gradual change in high-pressure level dependent on the motor rotational speed. 

Only one from four production six-hole injectors was installed onto the CRIS. The injector 
was set up horizontally into a suction duct to remove residual diesel spray during the 
measurements. The fuel firom both the common rail spill valve and the injector spill valve was 
returned back into the fuel tank through a flat plate water cooler. 

To control the high-pressure level into the common rail through its spill valve, a pressure 
limit control was employed in the system. A TTL type 200 Hz 10 V 70 % duty cycle voltage 
signal was coded into an arbitrary waveform generator by using bench link based software. An 
electronic limit switch controlled the final setup of pressure limit. This electric signal was 
transmitted to a voltage-to-current converter that was constructed by employing an insulated gate 
bipolar transistor with an ultra fast soft recovery diode. 

The waveform generator output signal was connected to a gate pin of the transistor. The 
collimator-emitter pins were powered by a triple output DC regulated power supply, the same 
type of power supply used for the pressure limit switch. Therefore, the CRIS pressure level was 
set up in three stages. First, the low-pressure pump was set at 20 bar (290 psi) just using a 
hydraulic control valve. Second, using the motor rotation speed control, pressure was increased 
up to 100 bar (1450 psi). Finally, increasing the voltage through the gate of the transistor, 
pressure was set at the desired level between 1200 to 1900 bar depending on the multiple 
injection profile (the number and duration of injection shots). 

ROSA type EDU 

To build up a ROSA EDU channel, the following sub-system has been designed, constructed 
and utilized on a production Bosch CRIS applied to E-class European passenger cars. A 
commercially available inductance L/C meter with resolution down to nH was used to measure 
inductance of each injector installed onto the CRIS. A second fimction/arbitrary wave generator 
was incorporated into the system in order to code ROSA type special voltage time series and 
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afterwards to have an output that represents multiple injection signals. A 500 MHz IGsa/s 
oscilloscope was applied to verify the quality and actual time phase setups of the output control 
signal directed to the CRIS injectors. 

The entire multi-steps and muhi-loop ROSA design algorithm of this embodiment can be 
divided into three large stages: 

First. The procedure begins from measurements of electric properties of the injector such as 
inductance L and resistance R , to evaluate time (or frequency) response. That allows a 
calculation of energy transferred per each transient fraction of each injection event. Calculating a 
predetermined ratio of the energy transfer, e.g., the integral energy generated by ROSA over the 
integral energy that was designed for this specific injector solenoid reflected into current-time 
profile, it becomes possible to calculate parameters of the secondary coil (ROSA) which 

must generate a transient current for rapid operating of the valve. 

Second . In the next stage, one needs to construct a so-called "I-Function" current as a timely 
fractional series and determine a charging time interval that is applicable for rapid and stable 
control over the injector. An example of the I-Function shape is shown in Fig. 32. For internal 
combustion injectors with an electronically controlled hydraulic valve, at the valve opening 
stage the most critical part within given time interval is a fraction from the beginning of the 
injection profile to a phase where I-Function current reaches maximum because instantaneous 
velocity of the solenoid armature is proportional to instantaneous current u = i-yfl/m . 

On other hand, in the case of an air intake valve it is necessary to have the time series 
extended to the moment where the first derivative of current becomes almost zero. This is due to 
proportionality between instantaneous acceleration (force) and current 
derivative a = {dil dtY ^[Ljm . If ROSA is desired as firmware, at this stage the algorithm 

switches to fabrication of the ROSA electric circuit and its tuning upon a specified injection 
mode. If ROSA must be implemented as a code source, the algorithm continues to the third 
stage. 

Third. The I-Function current time series must be fitted to a standard waveform function 
available in an arbitrary (ARB) wave generator. After fitting the derived I-Function to the 
waveform function algebraically, it is necessary to construct different transient phases of the 
injection cycle including individual injection shots and their |as-fractions. Finally, constructed 
current code is transferred into the given ARB-generator that next controls the injection profile. 

The shots' profiles must be constructed for each engine mapping point according to the 
engine speed-load and emission control. A full combination of the multiple injection profiles 
forms a library of the injection different waveforms (LIW). Afterwards, the entire LIW must be 
transferred into an electronic injection-driving unit (EDU), which communicates with the main 
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vehicle electronic control unit (ECU). Depending on the driving conditions, the ECU calls either 
OEM's or LIW's code related to the particular injection situation. 

ROSA Bench Model 

It is necessary to know the exact operation data of a production injection system, for 
instance, injector current/voltage trace applied on its actuator. In the Bosch CRIS injector, the 
solenoid triggers a ball type valve. At the stage of its pulling in (energized solenoid) the bleed 
orifice is opened and pressure difference between the feed passage to the nozzle and valve 
control chamber causes upward lift of the nozzle needle sequentially resulting in the injection 
event. 

Typical current trace applied to the Bosch CRIS injector is illustrated in Fig. 33. The 
energizing time of this solenoid varies from 1 to 2 ms with a peak puUing-in current of 1 8A and 
holding current of 12A. The rise time and fall time are varied from 80 to 100 \xs. During the 
holding stage current oscillates with amplitude 0.57A and periodicity 0.1-0.2 ms. 

The power E = A [lI^)/a t fluxed into the primary solenoid during energized state is 
calculated using measured inductance L , pulling-in peak /^^^^ and holding /^^^^ current, time 
response and holding duration respectively A/ to peak and holding stages. Ep^^j^ varies from 64.8 
to 72. 9 W and E^^^^ = 4.7-6.1 W for various injectors. These power (energetic) values are limited 
by construction of the coil, i.e., its inductance L and currents 1^^^^ , /^^^^ upon dynamic time 

response. To make the solenoid function very rapidly it is necessary to have an increased energy 
that will be released in a very short time. 

The distance between the high-pressure injector inlet to its nozzle is about 0.1 1 m. The 
sound speed under common rail of 1600 bar is -1700 m/s, so the time of pressure propagation is 
about 65 fas. That implies a magnitude of time fraction that must be comparable with minimal 
rise/fall time of the actuator resulting in high cycle-to-cycle stability (repeatability) of the 
multiple injection profile. 

The secondary coil does produce a quick power release on the primary coil to facilitate both 
rising and falling transitions. In the right gray part of the table the first input is power ratio 
between of the injector coil and £^^^,^2 of ROSA coil Ep^^^2= P^peaki » where factor F is 

varied between 1.5 to 4.0 depending on the actuator type and its application. In this particular 
case, it is maximized because for multiple injection with a fine inductance (high response time) 
the effect of rapidness is associated with high power ratio F = 4.0 . That permits the calculation 
of inductance of the ROSA coil f {E^^^f^^, T^^^^^^ , I peaki) , 

Conversely, the ROSA coil has a slower time response T^^aki^ ^^peaki ' where 2.0 < A: < 5.0 . 
Once again, because multiple injection requires very quick response over both injection shot and 
dwell interval between these shots, factor k = 2.0 is minimized. That results in resistance value 
/?2= ^2 ^'^peaki * Now, havlng frequency responses of both injector and ROSA coils, one can 
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construct the I-Function current (as discussed in detail elsewhere in the present application). 

The I-Function current trace and its first derivative are shown in Fig. 32. Because Rl L data 

are of the magnitude order of kHz, the time scale is scaled out to ms. The maximum current peak 

corresponds to 0.047 ms which relates to the maximum velocity of the primary solenoid 
armature. That time duration is a time /^^^^ that should be given for the ROSA coil for its 

charging before the energy is transferred into the primary injector coil. 

Waveform generator hardware can reproduce a variety of the current traces called standard 
waveforms as well as their different combinations. That moves the algorithm to the next step, 
which is a translation of the I-Function current into available standard functions and the time 
phases into a number of points within the injection cycle. For instance, in the software used in 
this ROSA development, one cycle is equal to 16000 points (pts). For the rise and fall I-Function 
current most fitting shapes are rise and fall. In normalized form, the voltage amplitude V is 
equal I. So, a matching factor should be derived from the comparison of I-and ARB functions at 
rise and fall fractions. Each injection shot was divided into 3 main sub-phases: rise, holding and 
fall transitions. They were translated into absolute and arbitrary coordinates of time and voltage 
amplitude. 

Fig. 34 demonstrates an example of the output signal for a six-shot muhiple injection at 
engine speed of 3600 RPM, the cycle duration is 360 cam [deg]. Here, the beginning of each 
cycle is referenced by a stroboscope second channel signal. The "Main 1" 600 |is shot is set up at 
180° (top dead center - TDC). Before TDC there are the "Pilot" 400 ^is and "Pre-M" 400 
shots, i.e., during compression stroke. 

The dwell interval "Dwell 1" between "Pre-M" and "Main 1" is set up as 200 ^s, while the 
dwell interval "Dwell 2" between "Main 1" and "Main 2" is 500 ^is. The "Main 2", "After-M" 
and "Post" are during the combustion power stroke and exhaust stroke respectively, as was 
shown Fig. 28. 

Volt-to-Amp Converter 

Having voltage arbitrary waveform for multiple injection, one needs another voltage-to- 
current converter to power the injector. Therefore, the second injection control channel was 
constructed as shown in Fig. 29 and 30. A voltage type injection signal coded as described above 
and transmitted to an arbitrary waveform generator. This signal was transferred onto a voltage- 
to-current converter of the same type that was used for the pressure spill valve control. The 
signal from the waveform generator controlled the gate pin while the transistor collimator- 
emitter pins were powered by the DC regulated power supply. This entire algorithm can be 
written as a program that will produce coding of all phases and shapes to generate the necessary 
waveforms including I-Function rise and fall fractions and holding stage. In other words, a 
special library can be written in a compressed form for easy translation of this library into 
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hardware (EDU) for further "call" type functionality. On the other hand such a library provides a 
variety of physically manufactured secondary coil drivers for different automotive applications 
(injectors, valvetrains and other rapidly operating actuators). 

HIGH SPEED VISUALIZATION 

Three different high-speed techniques were used to visualize multiple injection dynamics. 
First, a film camera was used at a lower speed of 5,000 fps to document 5- and 6-shot multiple 
injections with a high spatial resolution and a high sensitivity. Evaluation of the liquid spray tip 
velocity resulted in a maximum speed of 250 m/s, which is below of the speed of sound -320 
m/s under normal ambient pressure and temperature in the laboratory room. However, it was 
obvious that during experimentation with diesel multiple injection the shock waves sound was 
clearly heard. 

Second, a very thorough study was carried out using a stroboscope "freezing" technique to 
learn what level of temporal resolution must be applied to see more transient fractions in the 
spray dynamics, especially at the beginning of each shot during multiple injections, as well as to 
estimate the delay between the electrical command signal generated from the waveform 
generator and the actual shot. This study has shown that a faction of a few 10 |xs equivalent to a 
high-speed visualization at a few 10,000 fps is essential to observe the spray dynamics. Delay 
time was estimated to be over 400 [is. 

Third, a high speed CCD camera with a speed up to 40,500 fps (24.69 |as /frame) was used 
to make numerous measurements in a wide range of setups of the injection repetition rate, 
number of shots, shot duration and dwell intervals at various spatial resolutions of the camera. 
Below, more details for each of these studies are described. 



Filming at 5,000 fps 

The setup for the filming is depicted in Fig. 35. The injector was mounted side-off through a 
glass wall of the protection box into the center of a 220-mm cylindrical black-wall duct in order 
to extract a residual mass of the spray into an exhaust hose connected to an external ventilation 
system. A US quarter of 24.76 mm was glued on the front black panel mounted just behind the 
injector nozzle tip in order to have a spatial scale on the observation disk. For illumination of the 
spray flow a laser channel was built up using a copper laser at 40 W output power. The pulse 
width was adjusted to 25 ns. An output beam of 25 mm was collimated by a 3320-mm plane- 
convex lens and redirected by a mirror to a 24-mm quartz rod in order to produce a laser sheet. 
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Inclination of the injection jets at 35^ to a vertical plane necessitated the use of such a thick laser 
sheet. A stroboscope was set up on a tripod to illuminate the beginning of each injection cycle. 
The injection ARB generator synchronized the cycle through a four-channel digital delay/pulse 
generator, which was used to set up the strobe light at any fixed time phase, i.e., to "freeze" the 
spray dynamics at this particular phase with very high temporal resolution available down to a 
Pico-second. 

For preliminary filming of the spray a high speed camera with an electronic control system 
was used. The camera was mounted on a tripod in the front position normal to the laser sheet at a 
distance of 300 mm and connected to its power and control units. A synchronization signal from 
the camera was fed back to the laser controller. At a camera speed of 5,000 fps, the acceleration 
time was 0.90 s from total filming time of 3.60 s for standard film length of 122 m. A high 
sensitivity film of 400 asa was used because the duration of the laser pulse was only 25 ns per 
each 200 ^is frame. 

Two films were made. The first one was filmed for six shots per injection cycle at an engine 
speed of 1,200 RPM. The second was filmed for five shots per injection cycle at an engine speed 
of 2,400 RPM. An example of visualization of 400 \xs Pre-Main (top raw), 600 ^is Main 1 
(middle raw) and 500 jis Main 2 (bottom raw) shots are illustrated in Fig, 36. An insufficiency of 
temporal resolution was observed due to the fact that the estimated spray tip velocity was less 
than sound speed. For example, the frame on top left shows a time phase of the beginning of 
Pre-Main shot. The length of each jet at this particular moment is twice the size of the reference 
coin, i.e., 49.52 mm. The frame duration is 200 |as. Therefore the estimated velocity is about 
247.6 m/s, below the speed of sound of 320 m/s. This fact contradicts what was heard (a 
supersonic sound) during run of the injection. 

Stroboscope "Freezing" Technique 

Afterwards, a special study was conducted and focused on the minimum temporal 
resolution needed for the measurements. The stroboscope light with a pulse width of 176 ^s and 
247 [IS at a repetition rate of 30 and 10 Hz, respectively, was gradually shifted along the cycle 
time phase. The delay generator was used to increment the shift at 100, 10 and 1 |is of time. In 
other words, a simulation of high-speed visualization was an equivalent to 10,000 and 100,000 
and 1,000,000 fps. The second increment was the most balanced in terms of the time 
consumption and resolution high enough to resolve the spray dynamics. 

Measurement of the jet length at the start of injection has shown that the spray tip velocity 
is over 360 m/s (supersonic). Increasing the number of shots per cycle from one to six, one can 
easily hear a very harmonic single tone sound becoming more and more husky under multiple 
injection runs because the shots are distributed in non-regular time intervals according to the 
multiple injection concept illustrated by Fig. 28. 
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The "voice" of multiple injection is very specific and can be recognized after getting some 
experience. At a repetition rate of 30 Hz, the frequencies of multiple harmonics are varied from 
30 to 1,600 Hz. Another important observation that came out of the stroboscope study is that at 
any frozen phase within a given injection shot one can see a very stable picture over many 
cycles. There is no oscillation of any part of the jets, neither in length nor shape nor density. 
That was the first clear indication that ROSA produces multiple injections with very high 
stability at all reasonable low, medium and high engine speed. 

Visualization at Higher Speed 

To monitor detailed diesel spray including the development of very initial transitions, a 
high-speed CCCD type digital video camera was adopted and used at various operational speed 
of 9,000 / 18,000 / 27,000 and 40,500 fps with spatial resolution of 256x128, 256x64, 256x64 
and 64x64 pixels per frame respective to the camera speeds. By increasing the speed, the study 
was mainly focused on initial single spray development in order to measure the spray tip 
velocity and delay of the injection shots relative to electronic signal setups as well as the exact 
dynamic duration of shots and dwell intervals between them, especially between Pre-Main 1 and 
from Main 1 to Main 2. The layout and photo view of the setup of the equipment is depicted in 
Fig. 37 and 38. The camera system includes (i) a compact camera mounted on a tripod with a 3D 
rotational traverse, (ii) a processor with a memory capacitor of 200 GB, and (iii) a lap top 
computer with a recording and post-processing software. The processor was connected both to 
the PC through Ethernet card and a video monitor. A trigger-in remote control was used to start 
the recording process. 

A 5 W argon laser continuously emitted a beam of 3 mm (488 and 514 nm wavelengths), 
which was re-directed through a mirror to a fiised quartz rod of 3.86 mm. Because the laser 
beam was not specially conditioned (coUimated) the final laser sheet thickness was about 12 
mm. This thickness was less than the 21 mm needed to cover the whole spray field in the duct 
because the jets were inclined at 35 degrees from the cutting laser vertical plan. However, it was 
larger than the space maintained by the camera at its high operational speed. 

The camera was mounted on a tripod in front of the injector nozzle tip at a distance of 1 80 
mm and slightly rotated at 25^ to capture the first jet counter clockwise from the direction of the 
laser sheet entrance. Again, the stroboscope was used to flash the injection cycle start. Using a 
light bulb and setup of the processor in "live" regime, the camera was carefully focused on the 
injector tip in such a manner that the quarter coin, which referenced spatial scale, was also 
clearly seen during flashing the stroboscope and the stroboscope together with the laser sheet as 
shown on photo A and B in Fig. 38. 

During high-speed visualization the laser beam was set up at 80 % of its peak power of 5 W. 
Multiple injections simultaneously with stroboscope flashes were run and the recording process 
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was started by the trigger-in signal. More than 20 films were recorded for various engine speeds, 
number of shots, variety of injection mapping setups and dwell intervals between Pre-Main 1 
and Main 1 shots. 

Treatment Process 

All recorded high-speed films were processed as sequential time-series. Fig. 39 illustrates 
an example of such series. It comprises 9 frames filmed during the Pilot shot of the six-shot 
injection cycle. The camera speed was 1 8,000 fps and the engine speed was set up at 2,400 
RPM. Because a thin laser sheet was used due to the lack of energy at the high-speed 
visualization only a portion of flight trace associated with initial phases in the vicinity of the 
injection nozzle was recorded. As shown at the enlarged fi-ame, a dark population of pixels 
presented in all digital films characterized the liquid jet tip. 

Within all injection events, four stages could be observed. During the first, a liquid jet is 
developed with supersonic speed that will be discussed later on. During the second, at the 
moment of closing the injector valve, the spray flow is detached fi-om the injector nozzle but 
some portion of liquid jet is still taking place. During the third, only the spray field can be seen. 
During the fourth, the diesel spray that inclined fi*om the vertical plan is moved out of the laser 

sheet and only its residual part is traced in the vicinity of the injector nozzle. The stroboscope 
flash indicated the start of each injection cycle A^^, . This frame was set up as zero time, which 
was used for subtraction for each other sequential frames = Nj^^ - N^, . The absolute time was 
calculated as a product of fi*ame duration and sequential fi"ame / = * 7}^^^^ = N I Camera Speed - A 
length of liquid jet tip Lj^^ projected on the vertical plan was measured against the coin scale. A 

post-injection length of the visualized jet from the beginning of spray to the liquid population 
^post was also measured. This length was almost constant during a few fi-ames and later it was 

decreased due to movement of the spray out of the laser sheet. Such a procedure allows an 
estimate of the lowest magnitude of the projected jet speed Vj^^ - Lj^, /tj^^ . This velocity is 

reflected in all processed data. The inclination of jet at angle a implies that projected velocity is 
Ujei = ^jet /cos(a*') . Because a thin laser sheet was used, the real jet tip velocity might be 

slightly higher. However, the measurement of exact jet speed velocity was not the main 

objective of this study. At the first stage of data processing, the main objective was to measure 
actual duration of each shot upon the length I^^, fi'om the beginning of the injection event 

until the moment when the spray was detached and to estimate the velocity that was supposed to 
be supersonic. The length L^^^^ and time ^^,of post-injection spray were also measured, so 

^post = ^posi Impost ' Because this length represents only the visual part of residual spray, this 

velocity becomes zero and even negative, just to characterize a post-injection fraction of the 
injection event. 

An example of liquid jet dynamics for a six-shot injection under engine speed of 1,200 and 
at camera speed of 18,000 fps is depicted in Fig. 40. First, one can see that all shots have 
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supersonic velocity. The end of injection in the velocity diagram is characterized by the fall 
crossing the ZERO line and the oscillation parts in negative zone are related to post-injection 
dynamics of the spray. The actual dynamic dwell interval between Pre-Main and Mainl shots is 
517 |is, between Main 1 and Main 2 it was 763 |is while the electronic setups were 300 and 500 
|xs, respectively. In this particular case, delay of the shot phases with regard to the electronic 
signals was about 500 |is. These aspects, i.e., the dynamic shot duration and delay, will be 
discussed in detail in the next paragraph. 

At the second stage special efforts were focused on cycle-to-cycle variation, in other words 
to estimate at which time fraction the variation can be detected. That was possible due to 
recording multiple injection events at different camera speeds. To analyze cycle-to-cycle 
variability, each injection setup was recorded as a series of sequential cycles. An example of the 
treatment process for the six-shot injection cycle monitored at the camera speed of 40,500 fps is 
illustrated in Fig, 41, Here, only four first injection shots, namely Pilot, Pre-Main, Mainl and 
Main 2 are plotted as 7 frame series for each shot (horizontal raw) in three sequential cycle 
series (vertical columns). Because the duration of the frame is 25.69 \is, the total time scale for 
seven frames plotted in Fig. 41 is 172.84 |is. However, all injection event data were processed 
until the moment when the jet was detached from the injector nozzle, i.e., the real duration was 
longer than shown in this figure. The main objective of the treatment was to analyze actual 
timing of shots' duration and its time phasing within each given cycle. That allowed analysis of 
factors of stability and time/phase delay with regard to the electronic timing setup shown earlier 
in Fig. 34. From Fig. 41 one can see, at least qualitatively, a high repeatability of the injection 
events in sequential cycle-to-cycle series for each shot. It can also be seen that the most "weak" 
injection characterizes the Pilot shot. The most "dense" injection, as expected, is seen during 
Main 1 and Main 2 events. 



RESULTS AND DISCUSSIONS 

Common Observations 

Cycle-to-cycle analysis has shown that even at a camera speed of 27,000 fps (time 
resolution of 37.04 |as) there is no cyclic variability in all physical data processed and analyzed. 
That is why for all further illustrations obtained at the highest camera speed of 40,500 fps data 
will be discussed. All data processed for each cycle were put into the cycle summary as shown in 
Fig, 42. On the left side of this table are data related to the electronic signals came out from the 
wave generator. On the right side are data obtained from the high-speed visualization record. 
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From this particular example one can conclude the following: 

1) The flow dynamic duration of each shot is shorter than was in the waveform setup. Duration 
of the Pilot, Pre-Main, After-Main and Post was equally setup to 400 |as, however, in real 
dynamics they have different duration varied from 173 p.s up to 222 ^s. The ARB duration of 
Main 1 and Main2 shots were 600 and 500 |is, respectively. During multiple injection they were 
shorten to 272 and 346 ^is. 

2) Controversially, the critical dwell intervals Pre-Main to Main 1 and Main 1 to Main 2 were 
increased from 200 to 518 |as (dwell 1) and from 500 to 691 \xs (dwell 2), respectively. 

3) All phases are shifted to about 400 |as. This delay is directly associated with the pressure wave 
propagation time in the common rail. Its equal to a fraction of the CRIS double length over 
sound speed of compressible diesel fiiel under such high injection pressure (over 1,400 bar). 

4) In terms of cam angle positioning at this high engine speed regime 3,600 RPM, there is quite 
small phase fraction well controlled during multiple injection. For instance, three injection 
events namely Pre-Main, Main 1 and Main 2 are laid within 21.9 degree while total these three 
shorts duration is 2.1 |is. 

Further studies were focused on three physical parameters important to characterize stability 
or controllability of the ROSA multiple injection: (i) the injection shots duration, (ii) the stable 
phasing of injection shots and (iii) the delay between the dynamic injection events and the ARB 
setups produced by the injection generator. All these data will be presented in absolute time 
scale and cam phases within cycle of 360*^. To make such analysis, all high-speed data filmed at 
40,500 ^s for 6-shot injection cycle at engine speed of 1,200 / 2,400 and 3,600 RPM were 
sorted per each three cycles for each injection case. 

Analysis of Short Duration 

The shots duration and its standard deviation along with ARB shots duration setups are 
shown in Fig. 43. Looking at this parameter in absolute time scale (2 top plots) and in camshaft 
angular position (2 bottom plots), one can conclude that: 

1) The higher engine speed, the longer injection duration actually generated from the injector. At 
higher engine speed the pressure, dropped during previous shot, has higher repetition rate to be 
recovered. 
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2) The shortest duration is deah with Pilot, Pre-Main and Post injection shots, 1 15, 178 and 140 [is 
in average at engine speed of 3,600 RPM, respectively. The longest shot duration is observed 
always at Main 2 event being 337 |xs at the same engine speed. 

3) High standard deviation of 38 |is belongs to Main 2, After-M and Post injection while almost 
ZERO deviation shots are Pilot and Mainl, especially at higher engine speed of 2,400 and 3,600 
RPM. 

4) Each duration in cam degree scale is well resolved between shots on specific engine speed. There 
is no instability regarding misfiring of the injector. The standard deviation for most cases lies within 

0. 2^ except Main2 and Post at high engine speed. 

Phasing of Injection Shots 

The phasing of shots and its standard deviation is summarized in Fig. 44. The top 2 plots are 
related to the absolute time scale, the bottom 2 graphs are presented in cam angular scale. Three 
points are important to outline here: 

1 . From the correlation diagram seen on third plot from the top, one can conclude that all injection 
events are delayed with regard to ARB waveform setups. Here, the vertical axis represents ARB 
setups; the horizontal is reveal to the actual phasing of the shots. Most long delay is suited for 
Main2 shot at high engine speed of 3,600 RPM. Instead of 183.96^ it becomes 196.09^ That is 
why for multiple injection control it will be necessary to make start of the injection events in 
advance to the phases that desired from the point of combustion control. To decrease phasing 
delay it is also possible to increase the CRIS pressure level. That would results in increased 
sonic pressure wave propagation, since shortening a time to recover a pressure loss from 
previous Pre-Main and Mainl shots. 

2. In general, actual phasing deviations are increasing with gradually increased engine speed. 
From the second (absolute time) and fourth (cam angular phase) plots all deviation data are 
clearly separated for the engine speed of 1,200 (red squares) to 2,400 (blue triangles) to 3,600 
(brawn cycles) RPM, respectively. 

3, Almost all shots are characterized by deviation of 14 jis, only at high engine speed the Main 1 , 
After-M and Post shots have deviation of 29, 25 and 29 ^is. In terms of cam degree, almost all 
deviations are laid within 0.2^ and maximum high engine speed phase fluctuation is about 0.3^. 
These data prove the high stability in the phasing of injection shots within the injection cycle. 

Critical Dwell Intervals 

The most critical control of dwell intervals between multiple injection events (shots) is dealt 
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with dwells between Pre-Main and Mainl (dwell-1), Mainl and Main2 (dwell-2). The are two 
physical phenomena that limit shortest dynamic dwell interval. The first is the time response 
constant of the injector solenoid. To get injection started, the injector solenoid needs a time 
^response = ^ / ^ determined by inductance and resistance of the coil, i.e., its design characteristics. 

For the Bosch CRIS injectors used in present study, this time is varied from 146 to 191 ^is. 

The second dwell shortest limit relates to a pressure recovery time needed after previous 

injection event and associated with double length of the common rail and sound speed (pressure 
wave propagation) t^^^^^^^ =2L/a.As discussed before, based on visualization measurements, 

this time is about 400 ^is. That is why the total transient dwell time t^^^^ > z^.^^^, + ^pressure is 
about 550 ^s. 

As example to such explanation, the processed data are reflected in Fig. 45. During the 
measurements, the dwell-1 and dwell-2 were setup by using ARB generator at 200 and 500 |is. 
The actual multi-injection dynamic dwells were measured by the high-speed camera with 
resolution of 24.69 ^s. As shown, the dwell-1 is varied from 494 to 543 at different engine 
speed with standard deviation between ZERO and 43 |as while dwell-2 is oscillated between 601 
and 716 iLis with deviation of 14 to 25 \xs. 

On two diagrams in the bottom part of Fig. 45, one can see that there is clear gradual 

separation of measured data depending on engine speed. The faster engine speed, the longer cam 

interval is needed for both dwell-1 and dwell-2. The longer absolute dwell time, the longer 

camshaft rotation will occur. In terms of camshaft degrees, the standard deviation is lower than 

0.3^ at the high engine speed of 3,600 RPM. 

To reduce pressure recovery time /^^^^^^^^ , one needs either to fabricate a new multi-sectional 

common rail with shorten length of each chamber connected individually to each injector (in- 
line common rail - inexpensive solution) or drastically increase of the pressure level, which 
ultimately results in increased density and since that the sound speed (high pressure pump - 
expensive solution). 



CONCLUSIONS AND FINAL REMARKS REGARDING THE PERFORMANCE 
EVALUATION OF A MULTI-BURST RAPIDLY OPERATING SECONDARY 
ACTUATOR ACCORDING TO AN EMBODIMENT OF THE PRESENT INVENTION 

In this study a ROSA-based diesel multiple injection test cell was constructed as a broad 
bench model that generated up to 6 shots with empirically proven high stability. This stable 
operation was evaluated over a wide range of the engine speeds varied from 1 ,200 to 3,600 
RPM. 

Up to six shots were produced with the shortest dwell setup between Pre-Main and Mainl 
of 200 |xs that was almost equal to the time response constant of the CRIS injector solenoid. 
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Moreover, the ROSA-based control system permits to generate more than 6 shots within 
injection cycle due to flexible setup of the current peaks released in ultra-shot time fraction. 

On the basis of high-speed visualization of the diesel multiple injection spray dynamics, the 
cycle-to-cycle timing variability, the stability in the shots duration is detected to be within 40 |as 
in absolute timing or 0.4^ in cam angle. The standard deviation of multi-shot phasing is not 
longer than 30 |is or 0.3^. The stability in cyclic variation of the shortest dwell intervals is also 
proven to be within 40 |as or 0.4° over entire range of the engine speed. Such high stability both 
in the timing of injection shots duration and dwell intervals and the phasing of injection events 
within sequential injection cycles is not currently demonstrated by using any other multiple 
injection techniques. A number of general technical conclusions and remarks came out from this 
study: 

1 . A third type of ROSA was constructed and applied to controlling of highly stable diesel 
multiple injection process. It was applied onto existing diesel injection system without any 
redesign of the original CRIS and injector unit. The ratio of the injector inductance to its 
resistance was very low; lower than for other type of the hydraulically/electronically controlled 
diesel injectors, the air intake valve and the gasoline injectors. That drafts the first principal 
conclusion that ROSA technique is applicable to numerous other devices where either rapidness 
(diesel multiple injection) or high cycle-to-cycle stability (gasoline injectors) or controllable 
almost zero sealing velocity (variable intake valves) are critical factors for the driving control. 

2. The performed timing limits are not associated with ROSA itself, but rather with a complexity 
of the high-pressure wave dynamics and multi-frequency hydraulics. During multiple injections 
with different dwell intervals between injection events a series of harmonics is presented in the 
common rail and injector oscillatory flows. 

The higher frequency of oscillation, the shorter length of pressure wave propagation occurs into 
pressure system. That necessitates a possible solution for decreasing delay by subdividing a 
high-pressure chamber, for instance common rail, into a series of short sections. 

3. The ROSA technique generates multiple injections with the stability of 40 -50 ^s, which is 
detectable at the high-speed of visualization at 40,500 fps. Even at the speed of 18,000 and 
27,000 fps, "instability" was not detectable. This level of stability is much higher that needed for 
injection and combustion control in automobile industry. For commercial implementation of 
ROSA, an electronic unit may be installed on the vehicle board to work in communication with 
its ECU. The code, obtained after tuning ROSA onto specified engine, may be either written into 
a remote chip (processor) or directly into OEM's ECU chip. Depending on the cost of the 
technology and engine type, the main advantage of ROSA is very stable phasing, dwelling and 
duration of multiple injection shots proved from cycle-to-cycle analysis. 
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II. Quantification of Instantaneous Diesel Flow Rates in Flow Generated By a Stable and 
Controllable Multiple Injection System 



INTRODUCTION 

The following now refers to a multiple injection technique according to an embodiment of 
the present invention that been applied to a common rail injection system (CRIS). This technique is 
based on a rapidly operating electromagnetic secondary actuator (ROSA) that generates transient 
current to control primary solenoid of the diesel injector with highly repeatable stability. Many 
advanced types of multiple injectors are designed by introducing a piezoelectric actuator. A control 
and test system was constructed to evaluate the ROSA multiple injection properties, particularly the 
instantaneous flow rates. The system has produced up to six shots per cycle under injection 
pressures of 120 to 180 MPa at repetition frequency from 10 to 30 Hz. An LDA-based system was 
applied to obtain centerline velocity into fuel feed pipe flow. The high-pressure flow passed 
through a specially fabricated transparent intersection. No artificially seeded particles were 
introduced into the flow. The data rate was high enough in order to accurately resolve cyclic-to- 
cyclic variation of injection shots. For each injection setup more than 1000 cycles were measured, 
sorted and processed to obtain angular resolved values of the flow rate, pressure gradient and 
integrated mass related to each individual injection event. The mass distribution per each shot can 
be accurately controlled by the ROSA system by means of the injection pressure, frequency and 
dwell/duration timing of the injection events. Applied instantaneous flow rate technique can be 
widely introduced for calibration and test of various high-pressure diesel multiple injection 
systems. 

Volumetric or mass flow rate measurements are among the most important measurements 
applied into many industries and engineering control systems. Particularly, in the field of fuel 
injection systems (FIS) employed to internal combustion engines, precise instantaneous fliel/air 
flow rate measurements provide control of equivalence ratio that determines following after 
combustion process. Variety of measurement techniques and apparatuses are used to obtain such 
information. For instance, a Bosch type fuel flow rate indicator, based on pressure wave 
propagation forward and back to a gauge sensor, is widely used for quantification of fuel amount 
generated by high-pressure gasoline and diesel FIS. Fewer studies are related to other types of fuel 
flow rate sensors, for example, based on a miniaturized hot wire anemometer, i.e., two thin film 
sensors to measure bi-directional flow, that was installed into the body of common rail injection 
nozzle. Now, the flow rate measurements become more valuable since introduction of various 
diesel multiple injection systems and technologies. The inventor has developed a unique method 
according to an embodiment of the present invention based on a laser Doppler anemometer (LDA) 
and applied it to a low-pressure (6 bar or - 100 psi) gasoline FIS, a gasoline direct injection (DI) 
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injection system which pressure was varied from 50 to 70 bar (~ 1,000 psi) using only a laminar 
flow solution due to a low oscillatory Reynolds number. 

The full solution including a part for turbulent transient injection flow has been described 
with regard to higher injection pressures up to 2000 bar 30,000 psi) and more that directly relates 
to diesel FIS. As it will be shown later, the full scope solution is also needed to measure complex 
flow dynamics in Dl-gasoline injection systems, for instance, equipped with swirl dual switch 
injector where ultra- fast spray dynamics characterizes by a superposition of jet and umbrella type 
substructures. 

There are two main objectives of this study. The first objective relates to instrumentation of 
an LDA flow rate meter (LDA FRM) and its application for various FISs such as a 4 bar gasoline, a 
100 bar servo-jet and a 1800 bar diesel. It will be shown that in gasoline application one needs to 
seed the fuel flow due to lack of oscillatory pressure level needed to generate naturally seeded 
scattering particles in the flow. For higher pressure, the system works without a need to seed the 
fuel flow. This phenomenon was firstly proved in normal-heptane FIS and now used in diesel# 2. 
The second object is continuation of the ROSA-controlled multiple injection system evaluation, 
which discussion was started above. Briefly, ROSA is a system that can be applied on any existing 
diesel injector equipped with a solenoid type actuator that controls injection active phase such as 
common rail (CR), electronic unit injector (EUI) or hydraulic electronic unit injector (HEUI). The 
same as in previous study, ROSA was employed to a CR based injection system (CRIS) and 
generated up to six injection events (shots) per each cycle. Integrated ROSA-CRIS system has 
demonstrated high stability and repeatability in multiple Injection patterns. Now, to quantify the 
fiiel amount injected per each individual injection event - active injection and passive injection, 
LDA FRM was newly constructed and applied to measure both cyclically averaged and time arrival 
time series to obtain the flow rate data. 

Details of the quantification are as described below: 

EXPERIMENTAL TECHNIQUES 

Flow Rate Measurement Method 

Initially, the method for measurement of instantaneous volumetric flow rate was developed 
for a laminar fast oscillating pipe flows. The analytical solution is based on three equations written 
with respect to a non-stationary flow from which three instantaneous values - velocity, pressure 
gradient and volumetric flow rate can be derived. The pressure gradient is superposed by a Fourier 
expansion to fit any arbitrary periodic flow: 
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dz 



(1), 



where conjugated CC. represent complex arguments of a given value. Taking into account linearity 
of the Navier- Stokes momentum equation on the pressure gradient term and using a superposition 
for each induced harmonics, the exact solution for velocity field can be found as 



+C.C. 



(2), 



where Taylor number ^^^^ R^^nlv defines partial velocity profile that responds to a particular 

oscillation " «" , ^ is inner pipe radius and *^ is kinematic viscosity. Normalized ratio of dynamic 

and viscous forces results in the viscous time constant ^ , being in present experiments a 

T 

few hundreds of ms. In other words, if harmonic period ^« = 2;r / « longer than ^ , the 
corresponding velocity profile will be fiilly developed as shown in Fig. 46, i.e., a parabolic one in 
laminar flow. Otherwise, it will be not fully developed and built up as a flat-flow with a strong 
shear stress at the pipe wall. An integration of velocity over a circular cross section yields the 
volumetric flow rate: 











4v tr 









(3). 



Now, for reconstruction of equations (1), (2) and (3) one needs to deduce harmonics 
<Po...pn> from a time series either of velocity or pressure gradient. In dependence on measurement 
point into pipe flow and temporal resolution essential to detect pipe flow transitions, different 
measurement techniques can be applied. Present technique is based on a centerline time-dependent 
velocity deduced from equation (2): 



-CC. 



(4). 



The velocity time series can be accurately obtained from LDA measurements that set up to 
a number of bins A^g^p within the injection cycle and transformed into Fourier expansion 



Uir^Oj) = ^/f{c,e"-+C.c) 



(5). 



That permits to compute unknown values of 
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(6). 



Capillary injection pipe flow includes short-time fractions when the injector opens and 
closes. Fast transient regime occurs at these moments and to reconstruct the transient flow 
dynamics a high temporal resolution is required. LDA-based flow rate metering technique meets 
this requirement. Basic limit of the method is dealt with the oscillation Reynolds number 
Re^ < 700 Stokes layer thickness ^ = yl2v/co ^ jhe injection systems related to gasoline 

(3-6 bar) and DI gasoline (50-70 bar) engines can be satisfactory measured using this laminar 
transient pipe flow model. 

In order to obtain accurate flow rate measurements in diesel FIS, more comprehensive 
solution of the Navier-Stokes equations for turbulent flow in a circular pipeline is required. The 
derivation of the turbulent flow rate solution has been fiiUy described. There, the continuity, z- and 
r-momentum, conservation equations, governing a 2D time-dependent, compressible, axially 
symmetric, elliptic, turbulent pipe flow with the only force due to pressure, are resolved with 
respect to Reynolds decomposition parts, the mean and fluctuation (pulsation) parts, of the axial 
u=^u + u' = u^ + u^^ + u' radial ^ = ^ ^ = + ^ velocity components, which are included to be 

r 

measured by LDA system with required temporal resolution, and diffusion -function potential 
^ = O + ^ jj^g present technique is related to the following four timing variables: 



• An injection cycle period 7"-. 10 ms. 

• A total injection duration ^ - 1 ms. 

• LDA cyclic phenomena measurement time span ~ where ^ 10"^, controlled by an 
electronic bin number generator, so 1 ^s, 

• A w V autocorrelation function delay ^ 1-100 |xs, i.e., it is over the measurement time span 
At 



For a short dynamic period ^ ^' , the integration of the given variable ^ matches to its 
fluctuation part of the total value ^^^^ . Wise versa, integration within a large time interval ^ T 
results in the mean part. The main criterion to determine clock-watch resolution is related to 
harmonic Stokes layer thickness <^ = V2K/«fi; = VvA//«;r <A^ ^^^^^ v djes^l kinematic 
viscosity (--2-4.5 mm^/s) and A is an optic fringe span (- 1-4 \xm) in the LDA beam intersection 
point. 

With respect to pressure gradient, three parts are also superposed, so that: 
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where ^« is the stationary portion, is the oscillating portion and ^"-^ is the fluctuation portion. In 
the full turbulent pipe flow transport equations, there are diffusion terms of the first, second, third 
and higher orders. However, for the high-pressure fuel injection pipe flow, the radial partial 
derivatives are as small as two or three orders of magnitude vs. the axial partial derivatives. 

Therefore, the first order of the pressure diffusion terms and has to be considered 
for the integration procedures. In other words, in order to obtain instantaneous volumetric flow rate 
over a pipe cross section in the direction of the pipe axis, it is necessary to integrate the « velocity 
component and turbulent velocity correlation ^^^^ projected on the same pipe axis as follows: 



4^ 



(8). 



This flow rate reflects an effective axial velocity composing four terms, i.e., a stationary 
part associated with , an oscillatory part associated with , a " -pulsation part associated with 

, and a - pulsation part, associated with . Expression for velocity measured on the 
centerline ^ ^ of the flow is: 



^ Av ^ 



(9). 



Accordingly, the experimentally measured centerline velocity time series may be expressed 
as the Fourier expansion: 

where switching in FFT summation is dependent on the following criteria: 



nefl.Ar,! if <5„ J— >10 A 



(11). 



Comparing equation (9) and (10) gives final expression for the pressure gradient series, 
which are needed to compute the instantaneous flow rate, expressed by the equation (8): 



(12). 



Therefore, two different FORTRAN-based programs according to the present invention 
were written with respect to laminar and turbulent oscillatory pipe flows. The output of this 
software permits obtain not only information about instantaneous volumetric or mass flow rates, 
but also pressure gradient and integrated (accumulated) fuel mass: 



which can be compared with a mass balance measurement to estimation accuracy of the LDA 
measurement (its optical alignment): 



LDA Flow Rate Stand and Test Flow Rigs 

The diesel flow rate test stand is schematically depicted in Fig 47. It consists from 4 
subsystems: (i) a testing fuel injection system (FIS), here specifically based on a BOSCH CRIS 
type, (ii) an electronic injection driving unit (EDU), here constructed as a ROSA-control system 
described in detail elsewhere in the present application, (iii) a commercially available laser Doppler 
anemometer (LDA) and (iv) the present inventor's software that reconstructs LDA output velocity 
data into instantaneous volumetric/mass flow rates. The high-pressure fuel delivery line is 
connected to a measurement intersection (MI) mounted between pressure source (pump or CR) and 
Injector. A capillary quartz pipe was installed into MI to have an access for the laser beams and the 
light scattered into the injection flow. 

Two different Mis were constructed for present injection tests. The design details to the first 
one are shown in Fig. 48. This MM worked under injection pressure up to 140 bar (-2,000 psi) 
and used in the present study for measuring flow rates generated by the gasoline and servo-jet type 
injectors. In this case the quartz pipe length was 300 mm, the factor of 100 times to its inner 
diameter of 3 mm that permitted to calibrate stand for both laminar and turbulent flows under 
transient injection as well as at steady state regimes, i.e., in very wide range of flow rates, very 
accurately due to fully developed flow profiles. Only two 0-ring sets into the MI-1 construction 
hermetically isolated the quartz pipe. The second intersection MI-2, photo of which is shown in 




Vlda P - 



itl mass balance 



(14). 



Fig. 49 (vertical steel MI-2 setup seen right to the pressure gauge), was designed for high pressures 
up to 2000 bar (--30,000 psi). The main part of MI-2 is a quartz pipe with inner diameter of 1 .90 
mm, outer diameter of 6.06 mm and length of 40.10 mm that was thermally pressed into a metal 
thick tube with outer diameter of 18,93 mm and length of 43.42 mm, designed and assembled 
according to the technique described earlier. 

Inner diameter of the cold steel tube before its thermal expansion at 600 C was 5.95 mm. 
So, after mounting the quartz piece inside of the heated tube and its slow gradual cooling, the 
quartz tube was strengthened due to radial strength from outer steel tube. That provided very good 
withstanding to diesel injection pressures. Afterwards, this pressed-fit unit was assembled into the 
housing using eight M8 screws and another larger size three well adjusted steel sections: in/outlet 
parts and supporting middle section with two large holes for penetration true of the laser beam and 
scattered light. All parts were precisely machined for matching each other in the length and contact 
disks diameter. MI-2 was used for the test of ROSA-CRIS muhiple injection system. To have a fine 
alignment, the MI was flexibly mounted onto a heavy metal frame with 3D alignment and 
adjustment mechanics. Ml-outlet was fiirther connected to the test injector. For instance, as shown 
in Fig. 49, the MI-2 housing with two 14 mm windows setup for laser beam penetration was 
installed between CRIS and injector fiiel inlet. MI was installed on the feed line in close vicinity to 
the injector. Particularly, in this case the total length between LDA measurement point, where two 
laser beams were intersecting into a vertical plane having the flow axis in, and the needle part of 
injector was 0.34 m. Taking into account that acoustic speed into highly pressurized fuel liquid is 
about 2000 m/s, the time delay in velocity series, proportional to the double length, is about 300 jas. 
This delay was validated during the measurements. 

A fully configured LDA system, depicted in Fig. 50, was used to measure centerline 
velocity into the injection flow. LDA itself composes an ion 120-mW laser, the transmitting and 
photo-receiving optics, a photo-detector unit, a 2-channel signal processor and a 3D traverse 
system, on which 310-mm transmitting and 400-mm receiving optics was mounted as illustrated in 
Fig. 49 and 50. 

The receiving optics was setup off-axis from the transmitting plane. Oflf-axis angle is 
always varied upon the fiiel and injection pressure. In the test of gasoline injection (law pressure of 
3-6 bar), when 5-^m aluminum oxide solid particles were seeded into the flow, any off-axis angle, 
even backscattering, was reliable to receive an LDA signal with high data rate. While diesel servo- 
jet diesel injection (medium pressure of 100 bar) was tested, the off-axis angle was set at 22^ after a 
number of alignment attempts. For ROSA-CRIS injection test (up to 2000 bar), it was found that 
39^ off-axis angle is the optimal for all measurement conditions. 

To monitor oscillatory injection flow, a cyclic phenomena type software was applied to sort 
and process LDA measurement data. To use it, an angular encoded startup signal was synchronized 
via a time delay generator by the same waveform generator, which controlled the injection duty 
cycle. The data rate was varied from 0.4 to 1 8 kHz that was enough to reconstruct multiple 
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injection cycle in all details of the magnitude and timely phased injection events. The LDA system 
measured velocity series in a reversible flow due to the electro-acoustic modulation (Bragg cells) in 
the transmition optics. Main parameters used for the measurements were: 



5. Phase averaging bins 360-3600 

Each centerline velocity time series were treated using the inventor's software. This 
program reconstructs the measurement data into instantaneous series of flow rate, pressure gradient 
and integrated (or accumulated) fuel mass within injection cycle. In order to determine whether 
laminar or turbulent flows are occurred during various injection runs, a variety of the flow rigs was 
studied: 

To simulate steady state flow, a water-filled vessel was elevated at different height. Under 
gravity force a seeded flow was streamed to a gasoline type injector that permitted to align the 
optical setup using max-velocity and min-rms criterion. 

A steady 10-bar pressurized water vessel, from which the fuel rail was connected to a 
gasoline injector. The measurements were obtained under pressure of 7.3 bar (-106 psi) at the 
injection frequency of 40 Hz. For this particular measurement the ROSA EDU was made as an 
electronic circuit sketched in Fig. 51. Only one control lag was used to facilitate opening of the 
injector valve. Two different ROSA secondary coil (SC) charging scenarios were applied as 
illustrated by Fig 52. Firstly, ROSA was charged from zero to 2000 microseconds and afterwards 
the primary solenoid (PS) in the injector was opened. The injection duration was the same for all 
measurements (15 ms). Secondly, the ROSA coil was charged fi*om zero to 2000 microseconds 
simultaneously with the injection signal applied to the primary coil. Injection duration was setup at 
3 and 5 ms, at each case a number of the instantaneous flow rate time series were measured. A 
combination of these two techniques results in phase-shifted or tuned charge scenario. 

A servo-jet type FIS was generated up to 100-bar pressure into delivering rail and up to 
1500-bar pressure in the injector accumulation branch. A stable LDA signal was obtained at the rail 
pressure over 40 bar. Non-seeded diesel # 2 fuel was. For measurements in the ROSA-CRIS 
multiple injection system, the injector, used in high-speed visualization, was mounted vertically 
onto the CRTS rail as shown in Fig. 47. Injector nozzle housing with diameter of 18.88 mm, was 
fixed inside of a metal tube connected in series with a pipe directed into a glass vessel to collect the 
injected fiiel settled on the mass balance. 

Calibration Procedure 



1 . Optical probe 

2. Fringe spacing 

3. Frequency shift 

4. Cyclic length 



size 77 X 77 X 945 \im 
3.15 jam 
40 MHz 



360^ 
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Simultaneously with LDA time series, an automated fuel mass data acquisition was run to 
obtain mean mass rate measurements accumulated into the vessel. The oscillating flows were 
measured in both laminar and turbulent areas. The results of comparison of the LDA and mass 
balance (MB) measurements in terms of mean velocity and mass rates are shown in Fig. 53. The 
split between laminar to turbulent zones lays at the mean velocity of 33 cm/s or the mean mass rate 
over 2 g/s. In laminar area the disagreement between LDA and MB is varied from -A to +2 %. In 
the turbulent zone it is shifted to -2 to 4 %. Integrated LDA system and software gives a good 
agreement, enough for calibration different FIS. The statistic correlation between LDA and MB 
measurements shown as the trend-lines in the figure indicates accuracy of 0.1 % for the mean flow 
rate in laminar flows and 0.7 % for the mean flow rate in turbulent flows. The total injection rates 
in ROSA-CRIS injection are more than 2 g/s, so only turbulent model is applicable to treat LDA 
velocity time series. Because different transient stages occurred during fuel injection as shown in 
Fig. 54, only linear "measured" part of the trace with the highest derivative was used for the final 
LDA-MB correlation. Data acquisition transient time was varied from a few seconds to a few tens 
of seconds dependent on the injection repetition rate, so more than a few hundred cycles were 
averaged during the mass balance measurement. 

In order to analyze and couple the fuel flow rates injected per each individual shot such as 
the Pilot, Pre-main, Mainl, Main2, After-Main and Post, the same multiple injection profiles those 
used before for the high speed diesel spray visualization were applied to the flow rate 
measurements. For each engine speed the original Bosch-type injection profile with duration of 2 
ms was also measured as a referenced fuel mass characterizing conventional CRIS operation. 

In Fig. 55 the data were measured at 30 Hz repetition rate for the referenced Bosch, ROSA 
single 600 |is shot and ROSA 6-shot injection cases. These are among most critical measurements 
because the high repetition frequency is associated with the high vibration of the fuel delivery line 
and pressure oscillation frequencies (30-1600 Hz). The disagreement between LDA and MB data is 
varied only in negative area from -11 to -4 %. 

In order to evaluate the fuel mass rates injected per each individual shot a mass extraction 
method was applied using only mass balance (MB) measurements. First, only one Mainl shot was 
generated by ROSA-CRIS system. The MB-time series was measured and the Mainl averaged 
injected mass ' was obtained. Second, the Pre-Main shot was added and a fiiel mass injected 
per two-shot injection cycles was measured. Since, the Pre-main injected mass was subtracted from 
current measurements " ^'^j " \ ^ This sequentially mass adding procedure was repeated 
until 6-shot injection profile was measured and last Post injection event was subtracted. Due to the 
problem of pressure recovery into CRIS, for the different engine speed the different pressures were 
generated: 1,600 bar at 1,200 rpm and 1,700 bar at 2,400 and 3,600 rpm. The Bosch type single- 
shot injection with duration of 1 ms was also measured as a reference. 
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RESULTS AND DISCUSSIONS 



Referring now to verification of injection system rapidness and its stability in timing, there 
is no guarantee regarding the timing response of the whole injector system as depicted in Fig. 56, 
even if the electric output signal from ROSA EDU indicates fast response. Direct application of 
ROSA in automotive field is related to the diesel and direct injection gasoline engines where a 
stratified charge of the fuel mixed with airflow determines the quality of combustion. 

According to the objectives, i.e., the LDA-based flow rate instrumentation and the ROSA- 
controUed multiple injection, the following results and discussions are separated into three sub- 
sections. The first two are related to the low- and mid-pressure FIS represented by the gasoline 
(ROSA-controlled) and servo-jet type injection systems to demonstrate capabilities of the 
instantaneous flow rate technique. The third is dealt with both objectives. 

Gasoline Type Low Pressure Injection 

The flow rate series obtained by using three different SC charge techniques reflected in Fig. 
52 are depicted in Fig. 57. All the data were measured under the same conditions: injection 
frequency 50 Hz, injection pressure 7.3 atm and SC charging time 2.0 ms. The right figure shows 
instantaneous volumetric flow rate series and the left plot depicts integrated (or accumulated) 
injected fuel mass. The first time series (black one) in both plots relates to simultaneously charging 
of the primary (injector) and secondary (ROSA) coils. The second line (red one) represents pre- 
charge scenario. The third curve (blue one) is the case when charging of SC (AC-wave form in Fig. 
52) has been started before the injection (CD-wave form in Fig. 52), however, at the moment of 1 .4 
ms when SC-charging was continued, the injection has been already run. So, the overlapping time 
was 0.6 ms. As one can see from instantaneous and integral time series, the fastest opening of the 
valve takes place under shifted (tuned) charge conditions. The slowest opening is associated with 
the pre-charge. This case also gives lowest level of flow amplitude meaning the lowest speed of the 
needle at the opening moment. A rapid response without any substantial phase delay is associated 
with the simultaneous charge of SC and PC. Essentially, the same flow amplitude characterizes 
both simultaneous charge and shifted charge. For diesel engines, where multiple injection events 
must be precisely phased and inject a larger amount of fuel, the shifted or "tuned" charge technique 
is mostly suitable. 

Details with respect to each charging scenario at the beginning phases (opening of the valve 
and startup of injection) are shown in Fig. 58. There are three plots of instantaneous volumetric 
flow rates at the top row and three plots of integrated (or accumulated) fuel masses at the bottom 
row. The first column reflects data obtained while SC was simultaneously charged with PC 
(injector), i.e. according to Fig. 51, i.e., ^ -timing was the same to ^-timing. The second column is 
related to measurements when SC was pre-charged before the injector PC (first was and 
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afterwards started ^ = ^in Fig. 51). The third column shows results when SC charging was 
shifted with respect to the injector PC operation, i.e., and intervals were overlapped. Under 
simultaneous charge, the longer the SC-charging time, the faster opening of the valve is observed in 
instantaneous series as the shift between different series towards the initial zero phase. The 
integrated mass series indicate increased speed of the valve that obviously seen through the slope 
S^^^S . In pre-charge case, increasing the charge time results in the same phase of the injection 
startup, but the amplitudes in the instantaneous series and the slopes in the integral mass series are 
gradually increasing that means increased injector valve speed. Both effects, the increased 
amplitude/slopes and rapidness occur under shifted charge shown in the third column of Fig. 58. 

Mid Pressure Injection (Servo-Jet / bkm) 

These measurements were objected to align hydraulic and optic systems in order to 
demonstrate LDA measurements without artificial seeding of the fiiel (diesel # 2). In Fig. 59 the 
time dependent centerline velocity and volumetric flow rate time series are plotted for two flows. 
The first (lower level) was obtained into seeded water flow while it was injected through a gasoline 
injector, p= 7 bar. The second (higher level) is related to injection generated by a servo-jet type 
system, p= 62 bar. 

The timing of injection cycle was the same: injection repetition rate of 11 Hz (equal to 
1,320 RMP) and duration of 15 ms. This simple comparison of different injection pressures shows 
that increased pressure is reflected by much more transient fuel flow before active injection phase 
(before the main rise slope), during injection (zigzag-type point in the rise indicating primary 
break-up into the ftiel spray, and rapid closing of injection - main fall slope), and after injection 
(post injection oscillations). The velocity and flow rates are increased in one order of magnitude. 
Next Fig. 60 is related to the servo-jet series of the pressure gradient and occurred into high- 
pressure fiiel upstream of the injector and integrated fijel mass injected per cycle. The fijel is 
flowing during entire cycle because it flows into return line while the injector triggering solenoid is 
de-energized. 

The injection transient dynamics can be characterized also in details related to specifically 
determined time/angular phases. As illustrated in Fig. 61, there are two parts of the interest. The 
first is when the injector valve is opening (4 points phased between 81^ and 94.5^) and the second is 
when the injector valve commanded to be closed (3 points phased between 130^ and 134.5^). On 
the bottom part of the picture one can see the dynamics of velocity profiles reconstructed for the 
same points. The opening process is performed by a series of a rapidly growing flat-type velocity 
shape in the central vicinity of the pipe flow and a shear stress at the pipe wall. Because the time of 
the transition is much shorter than viscous time constant, the velocity profile cannot reach a shape 
of the fijlly developed turbulent flow. The development process is continued, however the valve is 
closed. At that moment the velocity profile starts to be reversed at the wall and integration of the 
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profile over the pipe cross section in many cases might result a negative flow rate following by a 
series of the pressure post-injection oscillations. 

High Pressure Injection (Diesel) 

Estimated Multiple Injection Masses 

The fiiel masses measured for each injection event are illustrated in Fig. 62 as a function of 
dynamic camshaft cyclic phase obtained from high-speed visualization. A number of conclusions 
can be drawn down as the following. With increasing engine speed, the values of multiple and 
single Bosch-type injections are gradually increased. This fact is also true for the measurements at 
speed of 2,400 and 3,600 rpm where the average pressure in common rail was equal. The smallest 
fuel mass of 1.1 to 2.7 mg/cycle characterizes the Pilot shot. All sequential three shots, e.g., Pre- 
Main, Mainl and Main2, are increased with the engine speed, but at low speed the highest mass is 
related to Mainl . At higher engine speed Pre-Main becomes the dominant. Regarding two last 
shots, i.e., After-M and Post, at low engine speed the Pre-Main is higher than even Mainl and Post. 
Increasing speed, the Post injection is increased drastically. For illustration purposes, at the same 
cyclic phase as the Mainl, the integrated injected mass over entire 6-shot cycle and CRIS baseline 
single shot masses are also plotted. At low engine speed 1-ms referenced injection consumes 
almost two times more fuel (37.7 mg vs. 22.4 mg) than 6-shot multiple injection while the total 
actual duration of the last was 1 .8 ms. At medium and high engine speed the situation is reversed, 
i.e., the 6-shot injection results in larger mass than 1-ms single shot injection, mostly due to the 
increased mass of the Post. It means that the After-M and Post injection duration setups must be 
decreased from 400 ^s to 200 us that could result in a fiiel mass decreased in one order of 
magnitude. It is also important to outline that at the higher engine speed there is no need to have the 
After-M and Post injections. For instance, 4-shorts injection cycle consumes always fewer fuels 
than CRIS baseline injection cycle. The minimum measured value of injected mass is 1.2 mg, the 
maximum is 75.0 mg. 

The ROSA-based multiple injection control has very wide dynamic range, which is very 
important for practical application. Multiple injection dynamics is summarized in Fig. 63. On the 
top of plot, in order to have better readout resolution, the injected massed are plotted vs. angular 
phases coded as the electronic setups. As one can see, the increased engine speed increases the 
injection masses per shot per cycle. On the bottom part of the figure, the total 6- and 4-shot 
injection and the 1-ms CRIS baseline single shot injections are plotted as function of engine speed. 
At the higher engine speed not more than 4-shot injection is essentially needed for diesel 
combustion process. The fuel consumption ratio between 4-shot and single shot injections is 0.35, 
0.48 and 0.84 respectively to the engine speed of 1,200 / 2,400 and 3,600 rpm. 
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Frequency - Pressure Correlation 



The process of high-pressure oscillation in diesel FIS during multiple injection is very 
complex due to the essential setup of irregular dwell intervals between shots. According to the 
measurements, the shortest dwells were varied from 0.556 to 1.001 ms observed between Pre-Main 
and Mainl, Mainl and Main2, respectively. That results in a high frequency domain of 0.999 to 
1.799 kHz. Because other dwells between Pilot/Pre-Main, Main2/After-M, After-M/ Post, 
Post/Pilot are longer 1-10 ms), the low frequency domain varied from 0.021 to 0.253 kHz can be 
implied. It is different in one or two orders of magnitude with respect to the high frequency 
domain. Each harmonics reflects different time delay, pressure recovery time and reaction of CRIS 
to increased engine speed because each harmonic frequency is doubled or tripled by increasing 
injection repetition rate, but this multiplication factor is very different for the low and high 
frequency domains. High timing stability tested during high-speed visualization is due to very 
stable control of multiple injection in such comprehensive environment. 

Ratio of the injection duration of each shot ^ and dwell interval ^ suited before this shot 
plays a key role in control of stable injection. By relating each injection event to the factor of , 
the whole data are sorted into low and high frequency domains as shown in Fig. 64. The Mainl and 
Main2 high frequency injection events are varied in very small range because for a wider variation 
they will need higher pressure level to damp pressure distraction at these frequencies of - kHz. 
Reversibly, the low-frequency domain (Pilot, Pre-Main, After-M and Post) is very reactive to the 
change of any time scale, particularly dealt with engine speed at dwell interval of 3.498 ms (0.253 
kHz) related to Post injection at 3,600 RPM. It is also obvious that every shot has own resonator 
frequency indicated by a spike with increased injection fiiel mass at the medium engine speed. 

LDA Instantaneous Flow Rates 

Applied LDA system permits to measure velocity time series either upon time arrival of 
Doppier bursts (TA-series) or using cyclic phenomena by sorting data according to the cyclic phase 
within injection cycle (C-series). Obtaining TA-series is important to make a plan for the 
measurements under various injection timing and pressure conditions and to analyze cycle-to-cycle 
variability. To illustrate various measurement situations, three single injection TA-series are plotted 
in Fig. 65. The top of figure related to a low frequency injection 1.8 Hz, injection duration 10 ms, 
p= 1400 bar. In the mid, there is injection generated at frequency 3.2 Hz, duration 10 ms, p= 1800 
bar. At the bottom, the injection was produced at high frequency 110 Hz, duration of 3 ms, p= 1800 
bar. Along this order of diagrams, the data rate decreased from 3 kHz down to 51 Hz. That 
demonstrates that both, pressure and basic injection rate are very critical to have enough data to 
resolve injection transitions. 
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Pressure level gradually increases the data rate because increased intensity of the cavitation 
as expected. 

In next four figures Fig. 66 through Fig. 69 the measured data are presented as TA-series 
phased within the injection period (data rate --1-10 kHz). The following discussions are focused on 
four main output parameters produced by the processing code: (i) centerline velocity measured by 
LDA system, (ii) volumetric flow rate reconstructed through velocity and rms data using capillary 
pipe geometry and kinetic properties of the fuel, (iii) reconstructed pressure gradient and (iv) 
accumulated fuel mass. All data are correspondent to the injection cycle repetition rate is 10 Hz 
(1,200 RPM). In terms of camshaft, 1 ms is equal to 3.6° (100 ^is fraction is 0.36°). 

Fig. 66 illustrates injection dynamics generated by a 2-ms reference single injection. The 
start of injection (SOI) was setup at 1 80^, p= 1400 bar (--22,000 psi). One can see that before and 
after active injection, the entire dynamics is smooth enough. The injection shaped profile ends by a 
zigzag spike. The smoothness of this process is due to a low frequency of the pressure wave 
oscillation; basic oscillatory harmonic is 10 Hz. No other harmonics are occurred within the cycle 
and the time needed to recover pressure is long enough. Looking at accumulated fuel mass plot in 
Fig. 66, one can see that some of the fuel is flowing through the measurement intersection before 
and after active injection phase. Each injection event creates a local negative pressure gradient 
spike. After active injection, due to accumulated pressure in CR, the fuel flows towards injector 
through the feed pipe to balance the volume (mass) to be injected in next shot. This recovering 
balance will be discussed later with regard to 6-shot injection cycle. Its derivative (slope) increases 
with increased injection pressure, frequency and fuel mass. 

Fig. 67 represents dynamics for ROSA-controlled single injection, duration 600 ms, p= 
1600 bar. Here, it is possible to distinguish four different elements vs. previous lower pressure and 
long injection (2-ms single shot reference injection). 

Before and after injection there is relatively strong background oscillation that seemed 
initially like a measurement noise. However, comparing the accumulated mass series in Fig. 66 and 
Fig. 67, one can conclude that the higher pressure applied in this case causes the higher flow rate. 
The active injection duration itself characterizes by a cascaded profile meaning that the fuel spray 
is split into a number of the primary breakup like phases. Duration of the injection profile is 
obviously shorter than 2-ms injection profile shown in Fig. 66 as supposed. All values of the output 
parameters are increased due to increased pressure. 

In Fig. 68, ROSA-controlled six shots injection dynamics is presented by TA-series. The 
SOI setups for each injection event were 126^, 173^ 180^, 192^, 270^ and 315^ respectively to the 
Pilot, Pre-Main, Mainl, Main2, After-M and Post injection shots. According to the flow rate 
measurement, these phases are 126^ 175^, 182°, 186^ 270° and 315°. All events having long dwell 
interval before the shot are characterized by exact time/angular phase that was electronically setup; 
there is enough time to recover the pressure loss. Wise versa, in vicinity of 180° where three shots 
(Pre-M, Mainl and Main2) are setup closely (dwells 300 and 400 \xs), the phases are shifted 
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relatively to the initial SOI sets because pressure needs a time comparable with the delay constant 
(300 ^is). The sequential injection events can be well seen from the accumulated mass series 
represented by a cascade; the number of cascaded stages is equal to the number of injection shots. 

Fig. 69 shows details of all three injection series plotted together with a higher angular 
resolution. In velocities cycles, the peaks related to the referenced 2-ms single injection at 1400 bar 
has the same level that ROSA six-shot injection at 1600 bar, so the muhiple injection requires 
increased either the high-pressure level or dwell intervals for pressure recover. 

The peak flow rate per each shot is decreased during muhiple injection while the pressure 
increased up to 1600 bar vs. 2-ms single shot injection at 1400 bar. In the accumulated mass series, 
in multiple injection line one can see three flatted stages corresponding to the Pre-M, Mainl and 
Main2 events. 

For obtaining the fuel masses injected per each individual event during multiple injection 
shown in Fig. 68 the injection cycle was split into 1 1 intervals including 6 active and 5 passive 
injection intervals related to the injection and no-injection (recovering balance) stages. This 
instantaneous flow rate measurements were made with accuracy of -4.6 % according to eq. (14), 
i.e., mass measured by LDA system vs. direct mass balance rating. 

The results of integration are reflected in Fig. 70. Within accuracy of LDA measurements, 
the mass injected (38.17 mg) is almost equal to the mass (34.25 mg) that was delivered to the feed 
pipe (recovering balance). The smallest amount of fuel 4.18 mg was injected during Pilot shot, the 
largest 1 1 .65 mg was during Main2 shot. The cyclic resolution was setup at 360 bins per cycle. 
Increasing it to 3600 bins, the injection mass resolution can be about 1 |ig. ROSA control was set to 
resolve the wave form generation with resolution of 0.01 V, so increasing it to 0.001 V, the 
multiple injection control can resolve mass dosing at the level of 0.01 mg. Such level of control 
requires a high data rate over 10 kHz that can be technically reached at the injection pressure level 
> 1600 bar and injection frequency < 60 Hz (7,200 RPM). 

CONCLUSIONS REGARDING QUANTIFICATION OF INSTANTANEOUS DIESEL 
FLOW RATES IN FLOW GENERATED BY A STABLE AND CONTROLLABLE 
MULTIPLE INJECTION SYSTEM 

According to the two objectives stated above, the conclusions are also grouped into two parts: 
Instrumentation 

To test fuel dynamics generated by ROSA-controlled multiple injection system, a laser 
Doppler anemometer (LDA)-based system was constructed and applied to obtain instantaneous 
volumetric/mass flow rates measured in a CRIS-type diesel injection system and processed using 
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laminar and turbulent oscillatory pipe flow models. The high-pressure flow passed through a 
specially constructed transparent intersection in which press-fit steel-quartz tube cell was 
hermetically installed for introducing laser beams. No seeding particles were implemented for LDA 
measurements due to the nature of the high-pressure oscillatory pipe flow. High data rate permitted 
to resolve each injection event, i.e., its timing characteristics and masses distributed within 
injection cycle. Time arrival- and cyclic- type data were obtained and sorted upon the angular phase 
and processed to obtain time/angular resolved series of (i) flow rate, (ii) pressure gradient and (iii) 
integrated mass related to individual injections. This flow metering system was applied to a 
particular CR-type diesel injection system. But it is also applicable, for example, to any high- 
pressure FIS operating under injection pressure over 40 bar (600 psi): gasoline GDI- and diesel 
EUI- and HEUI-type systems. Such calibration stand can be used for the test, improvement, 
verification and certification of a variety of FIS components including injector itself. The technique 
provides wide dynamic range and high temporal resolution for flow rate measurements, including 
rapid transient reversible flow occurred during multiple injection cycle. 

ROSA Performance 

The mass rated measurements of individual fiiel masses injected during multiple injection 
controlled by ROSA-CRIS test system are shown promising results both in fuel dosing and 
injection control using low- and high-frequency domains associated with pressure wave 
propagation harmonics. 

The wide dynamic range (max-to-min) of the injected masses and well separated low and 
high frequency pressure oscillation domains provide a good validation for ROSA-type control in 
entire range of the engine speed, injection duration and setups of critical ultra-short dwells between 
injection events. ROSA injection control system produces highly stable phasing and duration of the 
multi-shot injection within 30 ^is as it was also detected by means of high-speed visualization of 
diesel sprays. The smallest mass injected is 4 mg, the largest is 18 mg. The mass distribution per 
each shot can be accurately controlled by ROSA system at the level as low as 0.5 mg by means of 
injection pressure, frequency and dwell/duration timing of the shots with the high measurable 
accuracy -0.01 mg. 

While a number of embodiments of the present invention have been described, it is 
understood that these embodiments are illustrative only, and not restrictive, and that many 
modifications may become apparent to those of ordinary skill in the art. For example, the code 
routines may be written in Fortran, a Fortran-like program, and/or any other program that will 
produce coding of all phases and shapes to generate special waveforms (including, for example, the 
I-Function rise and fall fraction). Further, a special library may be written (e.g., in compressed 
form) for easy translation library into hardware (e.g., an ECU) for further call type fiinctionality. 
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Further still, such a library may permit a variety of physically manufactured secondary coil drivers 
for different automotive applications (e.g., injectors, valvetrains and/or other rapidly operating 
actuators). 
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